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Abstract. The supporting worldsheet of a string, membrane, or other higher dimensional brane, 
is analysed in terms of its first, second, and third fundamental tensors, and its inner and outer 
curvature tensors. The dynamical equations governing the models appropriate for phenomena such 
as (superconducting) cosmic strings and cosmic domain walls are developed in a general framework 
(allowing for both electromagnetic and Kalb Ramond background coupling). It is shown how 
the surface stress momentum energy density tensor determines the propagation characteristics of 
small "wiggle" perturbations of the worldsheet. Attention is then focussed on special features of 
strings (using the transonic model with tension T inversely proportional to the energy density 
U as a particularly important example). A quadratic Hamilton- Jacobi formulation is shown to 
govern equilibium states and other conservative string configurations sharing a symmetry of the 
(gravitational, electromagnetic, and Kalb-Ramond) background, including stable ring states that 
may be cosmologically important. 



1 Worldsheet Curvature Analysis 

In preparation for the more specific study of strings in the last three sections of this course, the 
first three sections are intended as an introduction to the systematic study, in a classical rela- 
tivistic framework, of "branes" , meaning physical models in which the relevant fields are confined 
to supporting worldsheets of lower dimension than the background spacetime. While not entirely 
new[Q|g], this subject is still at a rather early stage of development (compared with the corre- 
sponding quantum theory || which has been stimulated by the rise of "superstring theory"), the 
main motivation for recent work[|] on classical relativistic brane theory being its application to 
vacuum defects produced by the Kibble mechanism J5|, particularly when of composite type as 
in the case of cosmic strings attached to external domain walls Q and of cosmic strings carrying 
internal currents [Q. 

Before discussing the dynamic laws governing the evolution of a brane worldsheet it is worth- 
while to devote this first section to a recapitulation of the essential differential geometric machinery^ 
needed for the analysis of a timelike worldsheet of dimension d say in a background space time 
manifold of dimension n. At this stage no restriction will be imposed on the curvature of the 
metric - which will as usual be represented with respect to local background coordinates x M (//= 
0, n-1) by its components g^ v - though it will be postulated to be flat, or at least stationary 



string and membrane theory has been hampered by a tradition of publishing results in untidy, 
highly gauge dependent, notation (one of the causes being the undue influence still exercised by 
Eisenhart's obsolete treatise "Riemannian Geometry" JlO[ ) . For the intermediate steps in partic- 
ular calculations it is of course frequently useful and often indispensible to introduce specifically 
adapted auxiliary structures, such as curvilinear worldsheet coordinates a 1 (i= 0, ... , d-1) and 
the associated bitensorial derivatives 

X •* da- ' { ' 

or specially adapted orthonormal frame vectors, consisting of an internal subset of vectors ( 
A= 0, ... , d-1) tangential to the worldsheet and an external subset of vectors A Y M ( X = 1, ... , 
n-d) orthogonal to the worldsheet, as characterised by 

1 a 11i bv = Vab j 'a^^xp = i ^x^yju = ^xy > (1-2) 

where rj AB is a fixed d-dimensional Minkowski metric and the Kronecker matrix S XY is a fixed 
(n-d)-dimensional Cartesion metric. Even in the most recent literature there are still (under 
Eisenhart's uninspiring influence) many examples of insufficient effort to sort out the messy clutter 
of indices of different kinds (Greek or Latin, early or late, small or capital) that arise in this way 
by grouping the various contributions into simple tensorially covariant combinations. Another 
inconvenient feature of many publications is that results have been left in a form that depends on 
some particular gauge choice (such as the conformal gauge for internal string coordinates) which 
obscures the relationship with other results concerning the same system but in a different gauge. 

The strategy adopted here[pd| aims at minimising such problems (they can never be entirely 
eliminated) by working as far as possible with a single kind of tensor index, which must of course be 
the one that is most fundamental, namely that of the background coordinates, x^. Thus, to avoid 
dependence on the internal frame index A (which is lowered and raised by contraction with the 
fixed d-dimcnsional Minkowski metric r/ AB and its inverse rj AB ) and on the external frame index 
X (which is lowered and raised by contraction with the fixed (n-d)-dimensional Cartesian metric 
8 XY and its inverse S XY ), the separate internal frame vectors t A M and external frame vectors A X M 
will as far as possible be eliminated in favour of the frame gauge independent combinations 

rf v = i A "i A u, = A X "A*„ , (1.3) 

of which the former, 77^, is what will be referred to as the (first) fundamental tensor of the metric, 
which acts as the (rank d) operator of tangential projection onto the world sheet, while the latter, 
_L M „, is the complementary (rank n-d) operator of projection orthogonal to the world sheet. 

The same principle applies to the avoidance of unnecessary involvement of the internal coordi- 
nate indices which are lowered and raised by contraction with the induced metric on the worldsheet 
as given by 
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and with its contravariant inverse K 13 . After being cast (by index raising if necessary) into its 
contravariant form, any internal coordinate tensor can be directly projected onto a corresponding 
background tensor in the manner exemplified by the intrinsic metric itself, which gives 

jf v = h^x^y^ , (1.5) 

thus providing an alternative (more direct) prescription for the fundamental tensor that was pre- 
viously introduced via the use of the internal frame in (1.3). This approach also provides a direct 
prescription for the orthogonal projector that was introduced via the use of an external frame in 
(1.3) but that is also obtainable immediately from (1.5) as 

-It = 9% - rf v • (1.6) 

In so far as we are concerned with tensor fields such as the frame vectors whose support is 
confined to the d-dimensional world sheet, the effect of Riemannian covariant differentation 
along an arbitrary directions on the background spacetime will not be well defined, only the 
corresponding tangentially projected differentiation operation 

%=V U ^, (1.7) 

being meaningful for them, as for instance in the case of a scalar field ip for which the tan- 
gentially projected gradient is given in terms of internal coordinate differentiation simply by 
V^V = h ij x^i ip :j . 

An irreducible basis for the various possible covariant derivatives of the frame vectors consists 
of the internal rotation pseudo-tensor p " and the external rotation (or "twist") pseudo-tensor 
as g iven by 

Pfp = rf„ i%% l a ° = ~P W V , w/p =± v „ \%%K a = V , (1.8) 

together with the second fundamental tensor K^ p as given by 

=± p a l a v % i A ° = -rf v \ x p %\ x <, . (1.9) 

The reason for qualifying the fields (1.8) as "pseudo" tensors is that although they are tensorial 
in the ordinary sense with respect to changes of the background coordinates x^ they are not 
geometrically well defined just by the geometry of the world sheet but are gauge dependent in the 
sense of being functions of the choice of the internal and external frames l a p and A X M . On the 
other hand, like the first fundamental tensor rf v as given by (1.5), the second fundamental tensor 
(1.9) is geometrically well defined in the sense of being frame gauge independent, as can be seen 
from the equivalent but more direct definition Q 

K^=n%%n p a . (1.10) 
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The gauge dependence of the rotation tensors p" and ui" p means that (unlike K^ p ) they 
can each be set to zero at any particular given point on the worldsheet by choice of the relevant 
frames in its vicinity. However the condition for it to be possible to set these pseudo-tensors to 
zero throughout an open neigbourhood is the vanishing of the curvatures of the corresponding 
frame bundles as characterised with respect to the respective invariance subgroups SO(l,d-l) and 
SO(n-d) into which the full Lorentz invariance group SO(l,n-l) is broken by the specification of 
the d-dimcnsional world sheet orientation. The inner curvature that needs to vanish for it to be 
possible for p^ p to be set to zero in an open neighbourhood is of Ricmannian type, is obtainable 
(by a calculation of the type originally developed by Cartan that was made familiar to physicists 
by Yang Mills theory) as[§ 

R kX » v = 2< CT 7 7 >' r [A V, t]/9 / T + 2p [lt ^p X]m> , (1.11) 

while the outer curvature that needs to vanish for it to be possible for the "twist" tensor w^ 1 ^ to 
be set to zero in an open neighbourhood is of a less familiar type that is given Q by 

n KX \ = 2 1.^ ± T fl 7f [A V K] ccv CT T + 2u [h ^u) X]w . (1.12) 

The frame gauge invariance of the expressions (1.11) and (1.12) is not immediately obvious, but 
will be made manifest in the the alternative expressions to be given below following a synopsis of 
the properties of the (by (1.10) manifestly) gauge invariant second fundamental tensor K pv p . 

An equation of the form (1.10) for K^p is of course meaningful not only for the fundamental 
projection tensor of a d-surface, but also for any (smooth) field of rank-d projection operators 
r\^ v as specified by a field of arbitrarily orientated d-surface elements. What distinguishes the 
integrable case, i.e. that in which the elements mesh together to form a well defined d-surface 
through the point under consideration, is the condition that the tensor defined by (1.10) should 
also satisfy the Weingarten identity 

K [fiv f = (1.13) 

(where the square brackets denote antisymmetrisation) , this symmetry property of the second 
fundamental tensor being derivable^ @] as a version of the well known Frobenius theorem. In 
addition to this non-trivial symmetry property, the second fundamental tensor is also obviously 
tangential on the first two indices and almost as obviously orthogonal on the last, i.e. 

±^K au P = K^ria" = . (1.14) 

The second fundamental tensor K^ v p has the property of fully determining the tangential deriva- 
tives of the first fundamental tensor r\ ^ v by the formula 

= 2K p{vp) (1.15) 
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(using round brackets to denote symmetrisation) and it can be seen to be characterisable by the 
condition that the orthogonal projection of the acceleration of any tangential vector field u^ 1 will 
be given by 

u^u v K^ p =± p ^ , it* 1 = u v V v u» . (1.16) 

It is very practical for a great many purposes to introduce the extrinsic curvature vector K^, 
defined as the trace of the second fundamental tensor, which is automatically orthogonal to the 
worldsheet, 

It is useful for many specific purposes to work this out in terms of the intrinsic metric hij and 
its determinant \h\. It suffices to use the simple expression V^i/? = h^x^^ip.j for the tangentially 
projected gradient of a scalar field ip on the worldsheet, but for a tensorial field (unless one is 
using Minkowski coordinates in a flat spacetime) there will also be contributions involving the 
background Ricmann Christoffcl connection 

= 9 V ° '{9<?(h,p) ~ \9y,p,o) ■ (1-18) 
The curvature vector is thus obtained in explicit detail as 

JC = %rT = -4= (VW\h ij x" ti ) d + Wx».x p .T» p . (1.19) 
V II I 

This last expression is technically useful for certain specific computational purposes, but it must be 
remarked that much of the literature on cosmic string dynamics has been made unnecessarily heavy 
to read by a tradition of working all the time with long strings of non tensorial terms such as those 
on the right of (1.19) rather than taking advantage of such more succinct tensorial expressions as 
the preceeding formula V^y^". As an alternative to the universally applicable tensorial approach 
advocated here, there is of course another more commonly used method of achieving succinctness 
in particular circumstances, which is to sacrifice gauge covariance by using specialised kinds of 
coordinate system. In particular for the case of a string, i.e. for a 2-dimensional worldsheet, it 
is standard practise to use conformal coordinates a° and a 1 so that the corresponding tangent 
vectors i M = x^ Q and = x 1 * satisfy the restrictions x^x'^ — 0, ± M ± M + x'^x'^ = 0, which implies 
yp^|| = x '^x\ L = -x^ so that (1.19) simply gives y/W\\ K " = x " u ~ x u + {x'^x'p - x»±P)T^ p . 

The physical specification of the extrinsic curvature vector (1.17) for a timelike d-surface in 
a dynamic theory provides what can be taken as the equations of extrinsic motion of the d- 
surface[4] [9], the simplest possibility being the "harmonic" condition = that is obtained 
(as will be shown in the following sections) from a surface measure variational principle such 
as that of the Dirac membrane model[Q, or of the Goto-Nambu string model(|] whose dynamic 
equations in a flat background are therefore expressible with respect to a standard conformal 
gauge in the familiar form x"^ — x^ = 0. There is a certain analogy between the Einstein vacuum 
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equations, which impose the vanishing of the trace TZ^u of the background spacetime curvature 
7Z\ti P v, and the Dirac-Gotu-Nambu equations, which impose the vanishing of the trace K v of the 
second fundamental tensor . Just as it is useful to separate out the Weyl tensorjl^], i.e. 

the trace free part of the Ricci background curvature which is the only part that remains when 
the Einstein vacuum equations are satisfied, so also analogously, it is useful to separate out the 
the trace free part of the second fundamental tensor, namely the extrinsic conformation tensor [pi, 
which is the only part that remains when equations of motion of the Dirac-Gotu-Nambu type are 
satisfied. Explicitly, the trace free extrinsic conformation tensor C Alt , p of a d-dimensional imbedding 
is defined[pj in terms of the corresponding first and second fundamental tensors n^ and K^p as 

V = K^p - -^p^KP , CV = • (1.20) 

Like the Weyl tensor Wj/» of the background metric (whose definition is given implicitly by 
(1.25) below) this conformation tensor has the noteworthy property of being invariant with respect 
to conformal modifications of the background metric: 

9^ ' * e 2a g^ , K^p i ^ K^p + Vfil/ If" V CT a , C^ p ^ C^" . (1.21) 

This formula is useful[^3) for calculations of the kind undertaken by Vilenkin|l4| in a standard 
Robertson- Walker type cosmological background, which can be obtained from a flat auxiliary 
spacetime metric by a conformal transformation for which e a is a time dependent Hubble expansion 
factor. 

As the higher order analogue of (1-10) we can go on to introduce the third fundamental tensor Q 

as 

S A/J! /= vW* ±P *VxK„ T a , (1.22) 

which by construction is obviously symmetric between the second and third indices and tangential 
on all the first three indices. In a spacetime background that is flat (or of constant curvature 
as is the case for the DeSitter universe model) this third fundamental tensor is fully symmetric 
over all the first three indices by what is interpretable as the generalised Codazzi identity which is 
expressible[|| in a background with arbitrary Riemann curvature 1Z\^ P a as 

Sa mi / = + f^r? V?" )ft CT A ± p (1.23) 

It is to be noted that a script symbol TZ is used here in order to distinguish the (n-dimensional) back- 
ground Riemann curvature tensor from the intrinsic curvature tensor (1.11) of the (d-dimensional) 
worldship to which the ordinary symbol R has already allocated. 

For many of the applications that will follow it will be sufficient just to treat the background 
spacetime as flat, i.e. to take TZ aT ^ a = 0. At this stage however, we shall allow for an unrestricted 
background curvature. For n> 2 this will be decomposible in terms of its trace free Weyl part 
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y^tiv p a 5which as remarked above is conformally invariant) and the corresponding background 
Ricci tensor and its scalar trace, 



n^ = n pil p u) n = n\, (1.24) 



in the form [h 2 



V = Wr" + lkK<\ - l^h^^A < (1-25) 
(in which the Weyl contribution can be non zero only for n> 4). In terms of the tangential 
projection of this background curvature, one can evaluate the corresponding internal curvature 
tensor (1.11) in the form 

= 2K p ^K v](JT +r ] ^r ] \ll^ T r ] p a if a , (1.26) 

which is the translation into the present scheme of what is well known in other schemes as the 
generalised Gauss identity. The much less well known analogue for the (identically trace free and 
conformally invariant) outer curvature (1-12) (for which the most historically appropriate name 
might be argued to be that of Schoutenflj^l) is given[|) in terms of the corresponding projection of 
the background Weyl tensor by the expression 

= 2C^ p C v]Ta + V > \W KX a r ±- p a ±- T „ ■ (1-27) 

It follows from this last identity that in a background that is flat or conformally flat (for which 
it is necessary, and for n> 4 sufficient, that the Weyl tensor should vanish) the vanishing of 
the extrinsic conformation tensor C^ p will be sufficient (independently of the behaviour of the 
extrinsic curvature vector K^) for vanishing of the outer curvature tensor fl^ p a , which is the 
condition for it to be possible to construct fields of vectors A^ orthogonal to the surface and such 
as to satisfy the generalised Fermi- Walker propagation condition to the effect that -L^V^Ap should 
vanish. It can also be shownQ (taking special trouble for the case d=3 ) that in a conformally flat 
background (of arbitrary dimension n) the vanishing of the conformation tensor C^ v p is always 
sufficient (though by no means necessary) for conformal flatness of the induced geometry in the 
imbedding. 

The application with which we shall mainly be concerned in the following work will be the case 
d=2 of a string. An orthonormal tangent frame will consist in this case just of a timelike unit 
vector, l p , and a spacelike unit vector, t^, whose exterior product vector is the frame independent 
antisymmetric unit surface element tensor 

E pv = 2 lo [A V ] = 2(-\h\y 1/2 x [p Q x y \ , (1.28) 

whose tangential gradient satisfies 

VxE" v = -2K Xp ^£^ p . (1.29) 
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(This is the special d=2 case of what is, as far as I am aware, the only wrongly printed formula in 
the more complete analysis S on which this section is based: the relevant general formula (B9) is 
valid as printed only for odd d, but needs insertion of a missing sign adjustment factor (— l) d_1 in 
order to hold for all d.) In this case the inner rotation pseudo tensor (1.8) is determined just by a 
corresponding rotation covector p^ according to the specification 

P A ^ = |f>A, px = p£ v £\. (1-30) 

This can be used to see from (1.11) that the Ricci scalar (1.24) of the 2-dimensional worldsheet will 
have the well known property of being a pure surface divergence, albeit of a frame gauge dependent 
quantity: 

R = %{£^p v ) ■ (1.31) 

In the specially important case of a string in ordinary 4-dimensional spacetime, i.e. when we 
have not only d=2 but also n=4, the antisymmetric background measure tensor £ X ^ u p can be used 
to determine a scalar (or more strictly, since its sign is orientation dependent, a pseudo scalar) 
magnitude for the outer curvature tensor (1.12) (despite the fact that its traces are identically 
zero) according to the specification 

fi= jfivp^" , (1-32) 

Under these circumstances one can also define a "twist" covector uj^ , that is the outer analogue of 
pfj,, according to the specification 

w„ = \u» x e Xy , pa £<>° . (1.33) 

This can be used to deduce from (1.12) that the outer curvature (pseudo) scalar ft of a string 
worldsheet in 4-dimensions has a divergence property of the same kind as that of its more widely 
known Ricci analogue (1.31), the corresponding formula being given by 

SI =%(£>»' w v ) . (1.34) 

It is to be remarked that for a compact spacelike 2-surface the integral of (1-29) gives the well 
known Gauss Bonnet invariant, but that the timelike string worldshccts under consideration here 
will not be characterised by any such global invariant since they will not be compact (being open 
in the time direction even for a loop that is closed in the spacial sense). The outer analogue of 
the Gauss Bonnet invariant that arises from (1.32) for a spacelike 2-surface has been discussed 
by Penrose and Rindlerp"5[ but again there is no corresponding global invariant in the necessarily 
non-compact timelike case of a string worldsheet. 
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2 Laws of Motion for a Regular Brane Complex 



The term p-brane has come into usc]3|] [[16| to describe a dynamic system localised on a timelike 
support surface of dimension d=p+l, imbedded in a spacetime background of dimension n>p. 
Thus at the low dimensional extreme one has the example of a zero - brane, meaning what is 
commonly referred to as a "point particle" , and of a 1-brane meaning what is commonly referred 
to as a "string". At the high dimensional extreme one has the "improper" case of an (n-1)- 
brane, meaning what is commonly referred to as a "medium" (as exemplified by a simple fluid), 
and of an (n-2)-brane, meaning what is commonly referred to as a "membrane" (from which the 
generic term "brane" is derived). A membrane (as exemplified by a cosmological domain wall) has 
the special feature of being supported by a hypersurface, and so being able to form a boundary 
between separate background space time regions; this means that a 2-brane has the status of being 
a membrane in ordinary 4-dimensional spacetime (with n = 4) but not in a higher dimensional 
(e.g. Kaluza Klein type) background. 

The purpose of the present section is to consider the dynamics not just of an individual brane 
but of a brane complex or "rigging model" jij such as is illustrated by the nautical archetype in 
which the wind - a 3-brane - acts on a boat's sail - a 2-brane - that is held in place by cords - 
1-branes - which meet at knots, shackles and pulley blocks that are macroscopically describable 
as point particles - i.e. 0-branes. In order for a a set of branes of diverse dimensions to qualify as 
a "geometrically regular" brane complex or "rigging system" it is required not only that the support 
surface of each (d-l)-brane should be a smoothly imbedded d-dimensional timelike hypersurface 
but also that its boundary, if any, should consist of a disjoint union of support surfaces of an 
attatched subset of lower dimensional branes of the complex. (For example in order qualify as part 
of a regular brane complex the edge of a boat's sail can not be allowed to flap freely but must be 
attatched to a hem cord belonging to the complex.) For the brane complex to qualify as regular 
in the strong dynamic sense that will be postulated in the present work, it is also required that a 
member p-brane can exert a direct force only on an an attached (p-l)-brane on its boundary or on 
an attached (p-l-l)-brane on whose boundary it is itself located, though it may be passively subject 
to forces exerted by a higher dimensional background field. For instance the Peccei-Quin axion 
model gives rise to field configurations representable as regular complexes of domain walls attached 
to strings Jl7|] fig] pS|] , and a bounded (topological or other) Higgs vortex defect terminated by a 
pair of pole defects |||| may be represented as a regular brane complex consisting of a 
finite cosmic string with a pair of point particles at its ends, in an approximation neglecting Higgs 
field radiation. (However allowance for radiation would require the use of an extended complex 
including the Higgs medium whose interaction with the string - and a fortiori with the terminating 
particles - would violate the regularity condition: the ensuing singularities in the back reaction 
would need to be treated by a renormalisation procedure of a kind(24j[^5| that is beyond the scope 
of the present article, but that is discussed by Paul Shellard in an accompanying article in this 
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volume and elsewhere[18].) 

The present section will be restricted to the case of a brane complex that is not only regular 
in the sense of the preceeding paragraph but that is also pure (or "fine") in the sense that the 
lengthscales characterising the internal structure of the (defect or other) localised phenomenon 
represented by the brane models are short compared with those characterising the macroscopic 
variations under consideration so that polarisation effects play no role. For instance in the case 
of a point particle, the restriction that it should be describable as a "pure" zero brane simply 
means that it can be represented as a simple monopole without any dipole or higher multipole 
effects. In the case of a cosmic string the use of a "pure" 1-brane description requires that the 
underlying vortex defect be sufficiently thin compared not only compared with its total length but 
also compared with the lengthscales characterising its curvature and the gradients of any currents 
it may be carrying. The effect of the simplest kind of curvature corrections beyond this "pure 
brane" limit will be discussed in Section 3, but in the rest of this article, as in the present section, 
it will be assumed that the ratio of microscopic to macroscopic lengthscales is sufficiently small for 
description in terms of "pure" p-branes to be adequate. 

The present section will not be concerned with the specific details of particular cases but with 
the generally valid laws that can be derived as Noether identities from the postulate that the model 
is governed by dynamical laws derivable from a variational principle specified in terms of an action 
function T. It is however to be emphasised that the validity at a macroscopic level of the laws 
given here is not restricted to cases represented by macroscopic models of the strictly conservative 
type directly governed by a macroscopic variational principle. The laws obtained here will also be 
applicable to classical models of dissipative type (e.g. allowing for resistivity to relative flow by 
internal currents) as necessary conditions for the existence of an underlying variational description 
of the microscopic (quantum) degrees of freedom that are allowed for merely as entropy in the 
macroscopically averaged classical description. 

In the case of a brane complex, the total action X will be given as a sum of distinct d-surfacc 
integrals respectively contributed by the various (d-l)-branes of the complex, of which each is 
supposed to have its own corresponding Lagrangian surface density scalar < d >£ say. Each supporting 
d-surface will be specified by a mapping a i— > x{a} giving the local background coordinates x M 
(/i=0, .... , n-1) as functions of local internal coordinates a 1 ( i=0, ... , d-1). The corresponding d- 
dimcnsional surface metric tensor ( d 'ft.y that is induced (in the manner described in the preceeding 
section) as the pull back of the n-dimensional background spacetime metric g^, will determine the 
natural surface measure, < d >(£>, in terms of which the total action will be expressible in the form 



d J 

As a formal artifice whose use is an unnecessary complication in ordinary dynamical calculations 
but that can be useful for purposes such as the calculation of radiation, the confined (d-surface 





(2.1) 
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supported) but locally regular Lagrangian scalar fields < d >£ can be replaced by corresponding un- 
confined, so no longer regular but distributional fields < d >£, in order to allow the the basic multidi- 
mensional action (2.1) to be represented as a single integral, 




over the n-dimensional background spacetime. In order to do this, it is evident that for each (d-1)- 
brane of the complex the required distributional action contribution < d '£ must be constructed in 
terms of the corresponding regular d-surface density scalar < d )£ according to the prescription that 
is expressible in standard Dirac notation as 

< d >£= \\g\\- 1/2 J^dS <- d >£5 n {x- x{a}} . (2.3) 

In the kind of model under consideration, each supporting d-surface is supposed to be endowed 
with its own independent internal field variables which are allowed to couple with each other and 
with their derivatives in the corresponding d-surface Lagrangian contribution < d >£, and which arc 
also allowed to couple into the Lagrangian contribution <- d - r >£ on any of its attached boundary 
(d-1) surfaces, though - in order not to violate the strong dynamic regularity condition - they 
are not allowed to couple into contributions of dimension (d-2) or lower. As well as involving its 
own d-brane surface fields and those of any (d+1) brane to whose boundary it may belong, each 
contribution < d >£ may also depend passively on the fields of a fixed higher dimensional background. 
Such fields will of course always include the background spacetime metric g^ itself. Apart from 
that, the most commonly relevant kind of backround field (the only one allowed for in the previous 
analysis [4]) is a Maxwellian gauge potential whose exterior derivative is the automatically 
"closed" electromagnetic field, 

F„„ = 2V [IX A V] , = 0. (2.4) 

Although many other possibilities can in principle be enviseaged, for the sake of simplicity the 
following analysis will not go beyond allowance for the only one that is important in applications 
to the kind of cosmic or superfluid defects that are the subject of discussion in the present volume, 
namely an antisymmetric Kalb-Ramond gauge field = —B vil whose exterior derivative is an 
automativally closed physical current 3-form, 

AW = 3V [AJ B„ p] , V {ll N vpa] = 0. (2.5) 

Just as a Maxwellian gauge transformation of the form i— > A^ + V M a for an arbitrary scalar a 
leaves the electromagnetic field (2.4) invariant, so analogously a Kalb-Ramond gauge transforma- 
tion i ► B^ u + 2V[ M x„] for an arbitrary covector \fi leaves the corresponding current 3-form 
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(5) invariant. In applications to ordinary 4-dimcnsional spacetime, the current 3-form will just be 
the dual Np Up — tp Vpa N' T of an ordinary current vector N 11 satisfying a conservation law of the 
usual type, V M iV = 0. Such a Kalb-Ramond representation can be used to provide an elegant 
variational formulation for ordinary perfect fluid theory |2(| and is particularly convenient for set- 
ting up "global" string models of vortices both in a simple cosmic axion or Higgs field||27||[ p8|p9|| 
and in a superfluid[pO| such as liquid Helium-4. 

In accordance with the preceeding considerations, the analysis that follows will be based on 
the postulate that the action is covariantly and gauge invariantly determined by specifying each 
scalar Lagrangian contribution ( d >£ as a function just of the background fields, Ap, Bp V and of 
course gp Ul and of any relevant internal fields (which in the simplest non-trivial case - exemplified 
by WittenjjJ type superconducting string models ||^] - consist just of a phase scalar ip). In 
accordance with the restriction that the branes be "pure" or "fine" in the sense explained above, 
it is postulated that polarisation effects are excluded by ruling out couplings involving gradients of 
the background fields. This means that the effect of making arbitrary infinitesimal "Lagrangian" 
variations d^Ap, d^Bp Vl d^gp V of the background fields will be to induce a corresponding variation 
dX of the action that simply has the form 

d1 = Y,\ <d) <^{^d,Ap + \ wW^B^ + \ WT"X<7^} , (2.6) 

provided either that that the relevant independent internal field components are fixed or else 
that the internal dynamic equations of motion are satisfied in accordance with the variational 
principle stipulating that variations of the relevant independent field variables should make no 
difference. For each d-brane of the complex, this partial differentiation formula implicitly specifies 
the corresponding electromagnetic surface current density vector (^J 11 , the surface vorticity flux 
bivector = -W' 1 , and the surface stress momentum energy density tensor wT^ v = 

WT U ^ , which are formally expressible more explicitly as 

_ S' d >C (5< d >£ — 6<- d lC — 

(d)jM _ _ ; (d)^" = 2 — , (.drp^ = 2 + < d >£ i^rf v , (2.7) 

5Ap ' 6Bp U ' 5gp V 

where CO^M" is the rank-d fundamental tensor of the d-dimensional imbedding, as defined in the 
manner described in the preceeding section. 

The condition that the action be gauge invariant means that if one simply sets d^Ap — V M a, 
d^Bp V — 2VruXvU d^gpv — 0, for arbitrarily chosen a and x M then dX should simply vanish, i.e. 

J2 I d^S { w7"V M a + id W^Vp X v } = . (2.8) 
d ^ 

In order for this to be able to hold for all possible fields a and x M it is evident that the surface 
current < d )J^ and the vorticity flux bivector WW**" must (as one would anyway expect from 
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the consideration that they depend just on the relevant internal d-surfacc fields) be purely d- 
surfacc tangential, i.e. their contractions with the relevant rank (n-d) orthogonal projector ( d )_L^= 
g^ v — < d >?7^ must vanish: 

( d)_L^ wjv = q ) (d)J _^ (d )^P = o _ (2.9) 

Hence, decomposing the full gradient operator as the sum of its tangentially projected part 
< d >V Jli = ^rf^y and of its orthogonally projected part ^'-L^ V M , and noting that by (2.9) the 
latter will give no contribution, one sees that (2.8) will take the form 

f {d, <£ { (d % ( < d >J M a + wW^xJ) ~ a (d) % J" - Xv (d % (€ F } = , (2.10) 

in which first term of each integrand is a pure surface divergence. Such a divergence can be dealt 
with using Green's theorem, according to which, for any d-dimensional support surface ( d >«S of a 
(d-l)-brane, one has the identity 



< d >dS < d % < d >J M = j (d -^dS < d >J M (2.11) 

where the integral on the right is taken over the boundary (d-l)-surface of <9< d ><S of < d ><S, and < d 'A M 
is the (uniquely defined) unit tangent vector on the d-surface that is directed normally outwards 
at its (d-l)-dimensional boundary. Bearing in mind that a membrane support hypersurface can 
belong to the boundary of two distinct media, and that for d< n-3 a d-brane may belong to a 
common boundary joining three or more distinct (d+l)-branes of the complex under consideration, 
one sees that (2.10) is equivalent to the condition 

j^dS^a^%MJ' 1 - Y <d) V (d y'') + Xv( lp) % lp) W tn/ - Y <d)A M (d W^) | = , (2.12) 

p J ' d=p+l d=p+l 

where, for a particular p-dimensionally supported (p-l)-brane, the summation "over d=p+l" is to 
be understood as consisting of a contribution from each (p+l)-dimensionally supported p-brane 
attached to it, where for each such p-brane, (d) A^ denotes the (uniquely defined) unit tangent vector 
on its (p+l)-dimcnsional support surface that is directed normally towards the p-dimensional 
support surface of the boundary (p-l)-brane. The Maxwell gauge invariance requirement to the 
effect that (2.12) should hold for arbitrary a can be seen to entail an electromagnetic charge 
conservation law of the form 

<p%(p'J AI = <d) V (d) ^ • ( 2 - 13 ) 

d=p+l 

This can be seen from (2.11) to be be interpretable as meaning that the total charge flowing 
out of particular (d-l)-brane from its boundary is balanced by the total charge flowing into it 
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from any d-branes to which it may be attached. The analogous Kalb-Ramond gauge invariance 
requirement that (2.12) should also hold for arbitrary \v can be seen to entail a corresponding 
vorticity conservation law of the form 

<p%<p>VF^= ^ < i )A (1 (#''". (2.14) 

d=p+l 

A more sophisticated but less practical way of deriving the foregoing conservation laws would be 
to work not from the expression (1) in terms of ordinary surface integrals but instead to use the 
superficially simpler expression (2.2) in terms of distributions which leads to the replacement of 
(2.13) by the ultimately equivalent (more formally obvious but less directly meaningful) expression 

V M (5>>>)=0 (2.15) 

d 

involving the no longer regular but Dirac distributional current that is given in terms of the 

corresponding regular surface current < d >J M by 

< d >> = \\g\\- 1/2 J^dS wj^d n [x - x{a}] . (2.16) 

Similarly one can if one wishes rewrite the vorticity flux conservation law (14) in the distributional 
form 

V M (^J d W^)=0 (2.17) 

where the distributional vorticity flux ( d W^ u is given in terms of the corresponding regular surface 
flux wW^ by 

= |M|" 1/2 J id) dS id) W^5 n [x - x{a}} . (2.18) 

It is left as an entirely optional exercise for any readers who may be adept in distribution theory to 
show how the ordinary functional relationships (2.13) and (2.14) can be recovered by by integrating 
out the Dirac distributions in (2.15) and(2.17). 

The condition that the hypothetical variations introduced in (2.6) should be "Lagrangian" 
simply means that they are to be understood to be measured with respect to a reference system 
that is comoving with the various branes under consideration, so that their localisation with respect 
to it remains fixed. This condition is necessary for the variation to be meaningly definable at all 
for a field whose support is confined to a particular brane locus, but in the case of an unrestricted 
background field one can enviseage the alternative possibility of an "Eulerian" variation, meaning 
one defined with respect to a reference system that is fixed in advance, independently of the 
localisation of the brane complex, the standard example being that of a Minkowski reference 
system in the case of a background that is flat. In such a case the relation between the more 
generally meaningfull Lagrangian (comoving) variation, denoted by c? L , and the corresponding 
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Eulcrian (fixed point) variation denoted by d E say will be given by Lie differentiation with respect 
to the vector field £ M say that specifies the infinitesimal of the comoving reference system with 
respect to the fixed background, i.e. one has 

<k-d B = & (2.19) 

where the Lie differentiation operator ££ is given for the background fields under consideration 
here by 

i£A„ = CV V A^ + A v \7^ v , 

(£g^ v = 2V ( „&,) • (2.20) 

This brings us to the main point of this section which is the derivation of the dynamic equations 
governing the extrinsic motion of the brancs of the complex, which are obtained from the variational 
principle to the effect that the action X is left invariant not only by infinitesimal variations of the 
relevant independent intrinsic fields on the support suerfaces but also by infinitesimal displacements 
of the support surfaces themselves. Since the background fields A pi B pvi and g pv are to be 
considered as fixed, the relevant Eulcrian variations simply vanish, and so the resulting Lagrangian 
variations will be directly identifiable with the corresponding Lie derivatives - as given by (2.20) - 
with respect to the generating vector field £ M of the infinitesimal displacement under consideration. 
The variational principle governing the equations of extrinsic motion is thus obtained by setting 
to zero the result of substituting these Lie derivatives in place of the corresponding Lagrangian 
variations in the more general variation formula (2.6), which gives 

[ (d, <£ { (d) J M i£A„ + \ i'f" i£B^ + \ (d) T M!y i£g^} = . (2.21) 

d 

The requirement that this should hold for any choice of £ M evidently implies that the tangentiality 
conditions (2.9) for the surface fluxes wjv and ( d W^ must be supplemented by an analogous d- 
surface tangentiality condition for the surface stress momentum energy tensor wT^ v , which must 
satisfy 

( d)_L^ (djf»P = o . (2.22) 

(as again one would expect anyway from the consideration that it depends just on the relevant 
internal d-surface fields). This allows (2.20) to be written out in the form 

J2 / <d) ^{f (Fph <d) ^+ k N P^ (d>W ^ - (d % wT p - A p <d) ^ (d) ^ - B pv < d % < d W^) 
d •* 

+ ^V p ( y e{Ap (d) l^+B pu ^W^+ (#" p ))}=0, (2.23) 
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in which the final contribution is a pure surface divergence that can be dealt with using Green's 
theorem as before. Using the results (2.13) and (2.14) of the analysis of the consequences of gauge 
invariance and proceeding as in their derivation above, one sees that the condition for (2.23) to 
hold for an arbitrary field £ M is that, on each (p-l)-brane of the complex, the dynamical equations 

(p % <p)T% = Cp)/, , <p)/ p = <p)/ p + «/„ (2.24) 

should be satisfied, where the contact force density < p >/ p exerted on the p-surface by other members 
of the brane complex is expressible as 

(p) /p = J2 <d)A A< <d)T % . ( 2 - 25 ) 

d=p+l 

while the surface force density (p)/„ exerted by the background fields is given by 

(p, 7p = F m cp)7" + \N P11V . (2.26) 

As before the summation "over d=p+l" in (2.25) is to be understood as consisting of a contri- 
bution from each of the p-branes attached to the (p-1)- brane under consideration, where for 
each such attached p-brane, ( d 'A p denotes the (uniquely defined) unit tangent vector on its (p+1)- 
dimcnsional support surface that is directed normally towards the p-dimensional support surface 
of the boundary (p-l)-brane. The first of the background force contibutions in (2.26) is of course 
the Lorentz type force density resulting from the effect of the electromagnetic field on the surface 
current, while the other contribution in (2.26) is a Joukowsky type force density (of the kind re- 
sponsible for the lift on an aerofoil) resulting from the Magnus effect, which acts in the case of a 
"global" string |f27j ]28|| though not in the case of a string of the "local" type for which the relevant 
vorticity flux (pW" will be zero. As with the conservation laws (2.13) and (2.14), so also the 
explicit force density balance law expressed by (2.24) can alternatively be expressed in terms of the 
corresponding Dirac distributional stress momentum energy and background force density tensors, 
tdjj'A"' anc j wj^j which are given for each (d-l)-brane in terms of the corresponding regular surface 
stress momentum energy and background force density tensors ( d >T /iI/ and (d) / p by 

{d) T^ = \\g\\- 1/2 J^dS ^T^5 n {x - x{a}} (2.27) 

and 

«/„ = \\g\\- 1/2 J < d >dS Mj^lx - x{a}} . (2.28) 

The equivalent - more formally obvious but less explicitly meaningful - distributional versional 
version of the force balance law (2.24) takes the form 

^(5}^"")=//* (2.29) 

d 
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where the total Dirac distributional force density is given in terms of the electromagnetic current 
distributions (2.16) and the vorticity flux distributions (2.18) by 

h = f p» E (d)> + h N p^ Y. {A)W ^ • ( 2 - 3 °) 

d d 

It is again left as an optional exercise for readers who are adept in the use of Dirac distributions 
to show that the system (2.24), (2.25), (2.26) is obtainable from (2.29) and (2.30) by substituting 
(2.16), (2.18), (2.27), (2.28). 

As an immediate corollary of (2.24), it is to be noted that for any vector field k p that generates 
a continuous symmetry of the background spacetime metric, i.e. for any solution of the Killing 
equations 

V (A A) = , (2.31) 
one can construct a corresponding surface momentum or energy density current 

(pfpM = wr^k„ , (2.32) 

that will satisfy 

<p% <»>p> = <d>A M (d) ^ + (p) 7 M fcM • ( 2 - 33 ) 

d=p+l 

In typical applications for which the n-dimensional background spacetime can be taken to be flat 
there will be n independent translation Killing vectors which alone (without recourse to the further 
n(n-l)/2 rotation and boost Killing vectors of the Lorentz algebra) will provide a set of relations 
of the form (2.33) that together provide the same information as that in the full force balance 
equation (2.24) or (2.29). 

Rather than the distributional version (2.29), it is the explicit version (2.24) of the force balance 
law that is directly useful for calculating the dynamic evolution of the brane support surfaces. 
Since the relation (2.29) involves n independent components whereas the support surface involved 
is only p-dimensional, there is a certain redundancy, which results from the fact that if the virtual 
displacement field £ M is tangential to the surface in question it cannot affect the action. Thus if 
(p'i.^ £ y = 0, the condition (2.21) will be satisfied as a mere identity - provided of course that 
the field equations governing the internal fields of the system are satisfied. It follows that the 
non-redundent information governing the extrinsic motion of the p-dimensional support surface 
will be given just by the orthogonally projected part of (2.24). Integrating by parts, using the fact 
that, by (1.6) and (1.15), the surface gradient of the rank-(n-p) orthogonal projector <p'_L^ will be 
given in terms of the second fundamental tensor (v)K llv p of the p-surface by 

<p% = - MK^S - MK/ V , (2.34) 
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it can be seen that the extrinsic equations of motion obtained as the orthogonally projected part 
of (2.24) will finally be expressible by 

( P )T^ Wfi^/ = Cp)_L^ (p)/^ . (2.35) 

It is to be emphasised that the formal validity of the formula that has just been derived is 
not confined to the variational models on which the above derivation is based, but also extends to 
dissipative models (involving effects such as external drag by the background medium fl3|| [fl4| f33|| 
or mutual resistence between independent internal currents). The condition that even a non- 
conservative macroscopic model should be compatible with an underlying microscopic model of 
conservative type requires the existence (representing to averages of corresponding microscopic 
quantities) of appropriate stress momentum energy density and force density fields satisfying (2.35). 
This ubiquitously applicable formula is just the natural higher generalisation of "Newton's law" 
(equating the product of mass with acceleration to the applied force) in the case of a particle. 
The surface stress momentum energy tensor, MT^ U , generalises the mass, and the second funda- 
mental tensor, MK P , generalises the acceleration. This can be seen from the fact that in the 
1-dimcnsional case of a "pure" point particle (i.e. a monoplole) of mass to, the Lagrangian is given 
simply by <*>£ — —to, so the 1-dimensional energy tensor is obtained in terms of the unit tangent 
vector u M (u^u^ = —1) as ( 1 )T A " / = mu^u" . In this zero-brane case, the first fundamental tensor 
is simply given by Wr)'* 1 ' = —u p u v , and the second fundamental tensor is given in terms of the 
acceleration = u v V v u^ as WK^ = u tli u v ii p . Thus (2.35) can be seen to reduce in the case of a 
particle simply to the usual familiar form mu p =( 1 '_L P „ 

3 Perturbations and Curvature Effects beyond the Pure 
Brane Limit 

Two of the most useful formulae for the analysis of small perturbations of a string or higher brane 
worldsheet are the expressions for the infinitesimal Lagrangian (comoving) variation of the first 
and second fundamental tensors in terms of the corresponding comovings variation d^g pv of the 
metric (with respect to the comoving reference system) . For the first fundamental tensor one easily 
obtains 

d L *r = -rrir<k9pa , « = rr±° u d L9pa (3.1) 

and, by substituting this in the defining relation (1.10), the corresponding Lagrangian variation of 
the second fundamental tensor is obtainedj35| as 

d L iV =± P A nlrf v d L T a \ + (2 ±^ K v) Tp - K^frf^g^ , (3.2) 

where the Lagrangian variation of the connection (1.18) is given by the well known formula 

d,T a \ = g xp (V {a d,g T)p - \v p d^ T ) . (3.3) 
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Since we are concerned here only with cases for which the background is fixed in advance so that 
the Eulerian variation d E will vanish in (2.19), the Lagrangian variation of the metric will be 
given just by its Lie derivative with respect to the infinitesimal displacement vector field £ M that 
generates the displacement of the worldsheet under consideration, i.e. we shall simply have 

d L g ar = 2V (ff £ T) . (3.4) 

It then follows from (3.3) that the Lagrangian variation of the connection will be given by 

rf L r CT \ = v ia v T) e-n\ lTT)p e , (3.5) 

where H ur is the background Riemann curvature (which will be negligible in typical applications 
for which the lengthscales characterising the geometric features of interest will be small compared 
with those characterising any background spacetime curvature). The Lagrangian variation of the 
first fundamental tensor is thus finally obtained in the form 

c^rT = -2^V">£ CT , (3.6) 
while that of the second fundamental tensor is found to be given by 

(2 L% K v) » - gP T K^) (V.r + V T £ CT ) . (3.7) 

It is instructive to apply the forgoing formulae to the case of a free pure brane worldsheet, 
meaning one for which there is no external force contribution so that the equation of extinsic 
motion reduces to the form 

T^K^p = . (3.8) 

On varying the relation (3.8) using (3.7) in conjunction with the orthogonality property (2.22) and 
the unperturbed equation (3.8) itself, the equation governing the propagation of the infinitesimal 
displacement vector is obtained in the form 

L\ - n\ y ^) = -K^J dJP" . (3.9) 

The extrinsic perturbation equation (3.9) is by itself only part of the complete system of pertur- 
bation equations governing the evolution of the brane, the remaining equations of the system being 
those governing the evolution of whatever surface current |34| and other relevant internal fields on 
the supporting worldsheet may be relevant. The perturbations of such fields are involved in the 
source term on the right of (3.9), whose explicit evaluation depends on the specific form of the 
relevant currents or other internal fields. However it is not necessary to know the specific form of 
such internal fields for the purpose just of deriving the characteristic velocities of propagation of 
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the extrinsic propagations represented by the displacement vector so long as they contribute to 
the source term on the right of the linearised perturbation equation (3.9) only at first differential 
order, so that the characteristic velocities will be completely determined by the first term on the 
left of (3.9) which will be the only second differential order contribution. It is apparent from (3.9) 
that under these conditions the equation for the characteristic tangent covector \^ sa Y wu l be 
given independently of any details of the surface currents or other internal fields simply[Q by 

T^XnX* = • (3.10) 

(It can be seen that the unperturbed surface stress momentum energy density tensor plays the 
same role here as that of the unperturbed metric tensor g^ u in the analogous characteristic equation 
for the familiar case of a massless background spacetime field, as exemplified by electromagnetic 
or gravitational radiation.) 

The linerised perturbation equation (3.9) will of course completely determine the evolution of 
the displacements by itself if there are no internal fields, i.e. in the case of the Dirac Gotu Nambu 
model determined by a surface Lagrangian of the trivial constant form 

< d >£ = -m d (3.11) 

where m is a constant having the dimension of mass (which would be of the order of magnitude 
of the relevant Higgs mass scale in the case of a vacuum defect arising from the spontaneous 
symmetry breaking mechanism of the kind most commonly considered^)) and d is the dimension 
of the worldsheet (i.e. d=l for a simple point particle, d—2 for a string, and so on). In this rather 
degenerate special case the surface stress momentum energy density tensor determined according 
to the prescription (2.7) will be given simply by 

The unperturbed Dirac Gotu Nambu equation of motion is thus obtained in a form that is given - 
independently of the dimension d indicated by the prefix ( d ) which may therefore be dropped - by 
the well known harmonicity condition that is expressible as the vanishing of the curvature vector, 

K» = . (3.13) 

The corresponding perturbation equation is obtained from (3.9) in the form 

As well as the compact tensorial version |35| given here, the litterature includes other equivalent 
but formally more complicated expressions ]3(| |37[ (involving reference to internal coordinates or 
surface adapted frames of the kind discussed at the beginning of Section 1) that generalise earlier 
work restricted to the hypersurface supported ("wall" or membrane) special casepgfi p9| pfl| [pd| . 
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For most practical physical purposes the most useful generalisations of the Dirac-Goto-Nambu 
models governed by (3.12) and (3.13) are those of the very general category governed by (3.8) and 
(3.9) which allow for internal fields such as the currents that can be used, not only to represent 
the Witten type superconductivity effect in cosmic strings, but also also to represent the effect of 
ordinary elasticity in terrestrial applications, such as the strings of musical instruments which were 
already a subject of scientific investigation, albeit at an empirical rather than theoretical level, in 
the time of Pythagoras. However before proceeding to the discussion of such internal field effects 
in the following section, it is of interest to consider how the simple Dirac-Goto-Nambu model can 
be generalised in a way that goes beyond the "pure" brane description characterised in the free 
case by (3.8) and in the presence of an external force by (2.35). The distinguishing property of 
a "pure" brane model of the kind considered in the previous section is the condition that the 
action depends only on the undifferentiated background fields g^, A^, B^ Vl but not on their 
gradients. One of the most familiar kinds of example in which this condition fails to hold is that 
of an electromagnetically polarised medium, whose action^ depends not just on the gauge field 

but also directly on the associated field F^ u itself. However for lack of time and space we shall 
not consider such electromagnetic effects in the present section but will consider only the simplest 
category of "geodynamic" brane models, meaning those in which, as in the "pure" brane models 
of the Dirac-Goto-Nambu category, the action depends only on the imbedding geometry of the 
worldsheet and not on any other external internal fields. The simplest such extension of the "pure" 
Dirac-Goto-Nambu model (whose action is proportional just to the surface measure which depends 
only on g^ v but not its derivatives) is based on a Lagrangian consisting not only of a constant term 
but also of terms proportional to the two independent scalars that can be constructed as quadratic 
functions of the first derivatives of the metric, namely K^K^ and K tivp K llvp . The inclusion of 
such "stiffness" terms has been suggested by Polyakov and others [Q, one of the main reasons 
being allowance for the deviations from the "pure" Dirac-Goto-Nambu description of cosmological 



strmgM^MM or domain wall[^|^[|9j|5j|5l| defects that one would expect to arise if the 



curvature becomes too strong. The kind of Lagrangian constructed in this way, namely 

< d >£ = -m d + b WK P ^K p ~ c ^K^ p WK^P , (3.15) 



where m, b, and c are constants, has recently been the subject of several investigations |32| |53| [[34| |35| |56 
It is convenient to use the abbreviation 

WS^p = b WK P - c WK* vp (3.16) 

which enables the Lagrangian (3.15) to be expressed in condensed form as 

< d >£ = -m d + ^Spup ^K^ p . (3.17) 

The variation needed for evaluating the change in such an "impure" brane action will thereby be 
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obtainable from the formulae above in the corresponding form 

||/i||- 1 / 3 d L (ll^l| 1/a Z) = dJL+~£rr<k giu , (3.18) 

(again dropping the explicit reference to the brane dimension d) with 

rf L Z = (E X ^K X » - 2 E x »K x vp )d h9flv + (2 E» Xv - E^)V x d,g^ . (3.19) 

Although it is still possible to construct a formally symmetric stress momentum energy density 
tensor of the distributional type, the presence of the gradient term on the right of (3.19) will 
make it rather pathological, with not just a contribution proportional to a Dirac distribution as 
in the "pure" brane case described by (2.27) but also with a contribution proportional to the even 
more highly singular gradient of a Dirac distribution's). In order to be able to continue working 
in terms of strictly regular surface supported field, it is necessary Q to deal with the gradient 
dependence of the action in such an "impure" brane model by having recourse to the use of a total 
stress momentum energy tensor T M „ of the no longer no longer symmetric canonical type that can 
be read out from the variation formula 

|W|- 1/2 d L (IH| 1/2 £) =T» v %iC + 2EP° V„r~ %(2E/»V P C) (3.20) 

that is obtained after substitution of (3.4) in (3.18) and (3.19). The regular but non-symmetric 
canonical surface stress momentum energy density tensor obtained in this way]56| is given by 

T"„ = f» v + t% , (3.21) 

where T^ u is symmetric and purely tangential to the worldsheet (as a regular "geometric" stress 
momentum energy density tensor would be) with a merely algebraic dependence on the second 
fundamental tensor, being given by 

f \ = T„» = T\n\ = £r/t - 2 E x vp K/» , (3.22) 

while the remainder, which is of higher differential order, is expressible in terms of the third 
fundamental tensor (1.22) as 

t% = T\ L\= 2 c S Y„ ~ 2 b H" A \, , (3.23) 
which can conveniently be rewritten with the higher derivative contributions regrouped in the form 

#i - 2(c - b)S"„ + c rfprrWerx ±\ , (3.24) 

where 

E% = R\\ =± X V ^K X . (3.25) 
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It is to be noticed that the total canonical surface stress momentum energy density tensor obtained 
in this way is is still automatically tangential to the worldshheet on its first (though no longer on 
its second) index, i.e. 

± X „T* V = 0, (3.26) 

and that the higher derivative contribution proportional to the trace S M „ of the third fundamental 
tensor will drop out if the coefficients have the same value, c — b. 

The application of the variation principle to the effect that the surface integral of the variation 
(3.20) should vanish for any displacement £ M within a bounded neighbourhood can be seen to lead 
(via an application of Green's theorem as in the preceeding section) to dynamical equations |5t| of 
the form 

%T» u = 2 2 i rK IK ,» v . (3.27) 

As in the pure Dirac-Gotu-Nambu case discussed above, the foregoing system of equations is 
partially redundant: although it involves n distinct spacetime vectorial equations, only n-d of 
them are dynamically independent, namely those projected orthogonally to the d-dimensional 
worldsheet. The others are merely Noether identities which follow independently of the variation 
principle from the fact that a displacement that is purely tangential to the worldsheet merely 
maps it onto itself and thus cannot affect the action, as can be verified directly using the generalised 
Codazzi identity (1.23). 

The lack of symmetry of T% means that the construction of the corresponding momentum 
current vector, say, associated with a generic background spacetime Killing vector field will 
not be quite as simple as in the "pure" brane case for which an expression of the simple form (2.32) 
suffices. However using the well known fact that the Killing equation (2.31) entails the integrability 
condition 

V M V„fc" = K p „ pX k x , (3.28) 

together with the observation that the antisymmetric part of the canonical stress momentum energy 
density tensor (3.22) is given according to (3.23) just by 

T [H = _2V A £ A[ ^ ] , (3.29) 

it can be seen that, for any solution of (2.31), the ansatzj56| 

V p = T^k v + 2E"%V 1/ fc p (3.30) 

provides a surface current, 7"*, which satisfies the tangentiality condition 

±\ P"=0, (3.31) 

and for which the strict surface conservation law, 

%V» = , (3.30) 



23 



will hold whenever the equation of motion (3.25) is satisfied. It is to be remarked however that for 
a Killing vector of the irrotational kind for which V M fc M vanishes altogether the second term (in- 
terpretable as a surface spin density contribution) in (3.30) will not contribute, i.e. an expression 
of the simpler form (2.32) will suffice. This applies in particular to the case of an ordinary trans- 
lation generator in flat space, for which the corresponding conserved surface current will represent 
ordinary energy or linear momentum, whereas in the case of angular momentum the extra (spin 
density) term in (2.32) is indispensiblc. 

The preceeding formulae all include allowance for arbitrary background curvature, but, to 
obtain the analogue of the non redundant version (3.8) of the equations of motion in a reasonably 
simple form, the restriction that the background spacetime be flat, i.e. TV 1 vpa = 0, will now be 
imposed. This enables the required system of dynamical equations to be expressed in the form 

T^K^ p = 2(b - c) ±» p % ( lf v V*X) , (3.31) 

with the higher derivative terms grouped on the right hand side, which vanishes if b = c. 

In the particular case]5l|| of a membrane, meaning a brane supported by a hypersurface of 
dimension d=n-l, the second fundamental tensor and its trace will be given in terms of the unit 
normal A p (which in this case will be unique up to a choice of sign) by K^p — Kp V \ p and 
K p = KX p with K = K^t where K pu is the second fundamental form (whose sign depends on 
that orientation chosen for the normal). In an ordinary 4-dimensional spacetime background, this 
membrane case corresponds to d=3, for which the symmetric tangential part of the surface stress 
momentum energy density tensor will be therefore be obtainable from (3.22) in the form 

( 3 )jV = _ ( m 3 _ b < 3 >K 2 + c ^K p a ^K a p ) ~ 2b ^K^ + 2c ^K x>1 <&K X U . 

In the case of a string with d=2 one can use the fact that the trace free conformation tensor (1.20) 
will satisfy 2 (2) C A pp <2) C A / = < 2 >C k a p < 2 C kAp ^ 2) rjp U to obtain a corresponding formula (which holds 



regardless of the background spacetime dimension n) given 1 56 by 

(2 yp»» = _ m 2 (2)rj ^ + 2 ( c _ b ) mcv .vp (2)Rp _ 

4 Conservative String Models as examples of perfect Branes 

It is reasonable to postulate that a "weak" energy condition of the kind formulated and justified 



by Hawking and Ellis 1 57 should hold for any pure p-brane model as a condition for physical 
realism as a macroscopic description of a (p+l)-surface supported physical system at a classical 
level, meaning that the model's surface stress momentum energy density tensor (as introduced in 
section 2) should be such that the contraction T^ v fi^flv is non negative for any vector /3 M that 
is timelike. Furthermore the causality condition to the effect that there should be no timelike 
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characteristic covector (i.e. no superluminal propagation) can be seen from (3.10) to entail the 
further requirement (going marginally beyond the "weak" condition of Hawking and Ellis) that 
f3fj,/3 v should be strictly positive if f3^ is timelike. This leads to the formulation of what may 
be called the "minimal" energy condition for a pure p-brane which is expressible as 

<0 T^/y?„>0. (4.1) 

A (pure) p-brane model will consequently be characterised by a well defined surface energy 
density, U say, that is specifiable by an eigenvalue equation of the form 

= -17/9" (4.2) 

where the corresponding eigenvector /3 M is distinguished by the requirement that it be tangential 
and non-spacelike: 

= o , ^ ■ ( 4 - 3 ) 

It is apparent that the "minimal" energy condition (4.1) requires that the eigenvalue U should be 
strictly positive unless /3 M is null in which case it may vanish: 

U > , U = => = . (4.4) 

In the Dirac-Goto-Nambu model that is most familiar to present day cosmologists, the eigen- 
vector /J" is indeterminate and the energy density U is the same (in relativistic units such as are 
used here, with the speed of light set to unity) as the corresponding (surface) tension. However in 
general it is essential to distinguish the concept of energy density U from the concept of the tension 
(as used in physics since the formulation of Hooke's law at the time of Newton) from which term 
"tensor" is derived. The (surface) tension scalar, for which we shall use the traditional symbol T, 
is defincable generically in a manner consistent with traditional usage, for a "pure" p-brane (i.e. a 
(p+l)-dimensionally supported system) of the kind considered here by decomposing the trace of 
the surface stress momentum energy density tensor in the form 

T» fl = -U-pT . (4.5) 

Appart from the degenerate Dirac-Goto-Nambu case for which U and T are actually equal, the 
simplest possibility is that of a perfect p-braneQ, meaning one whose surface stress momentum 
energy density tensor is spatially (thus p-dimensionally) isotropic so that it is will be expressible, 
for a suitable choice of the normalisation of the eigenvector /?" in (4.2) by 

= - Trf v , (4.6) 

where the required normalisation is given by 

U — T = frp^ < . (4.7) 
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This category includes the case of the Dirac-Goto-Nambu model (3.12), which is obtained (with 
T = m p+1 ) by normalising the (in this case indeterminate) eigenvector to zero, i.e setting (3^ = 0. 
A more mundane example is provided by the familiar "improper" case p=n-l (where n is the 
background dimension) of an ordinary (relativistic) perfect fluid with pressure P = —T. Although 
the tension T is negative in the ordinary fluid case, it must be positive, as a condition for stability, 
for "proper" p-branes of lower dimension, p<n-l. Whereas the membrane models (with p=2 for 
n=4) that are appropriate for the description of nautical sails will not in general have the isotropic 
form (4.6), a familiar everyday example of a membrane that will automatically be "perfect" in 
this sense is provided by the case of an ordinary soap bubble whose boundary hypersurface will 
be characterised by < T « U. The kind of application on which we shall concentrate in the 
following sections is that of a string, as given by p=l, for which the "perfection" property (4.6) 
will always hold: isotropy cannot fail to apply in this 1-brane case because only a single space 
dimension is involved. 

Before specialising to the (automatically perfect) string case, it is to be remarked that except 
when the /3 M is indeterminate (as in the Dirac-Goto-Nambu case), or null (as can occur in special 
subspaces where the current becomes null in a string model of the kind ]3l|] p2| ]34|] appropriate 
for describing cosmic vortex defects of the "superconducting" kind proposed by WittenQ) this 
eigenvector will determine a corresponding unit eigenvector, 

«" = (U - T)- 1/2 (3» , - -1 , (4.8) 

specifying a naturally preferred rest frame, in terms of which (4.6) will be rewriteable in the 
standard form 

= {U- T)u»u v - Tif" . (4.9) 

With respect to the frame so defined, it can be seen that the velocity c E say of propagation of ex- 
trinsic perturbations (what are commonly referred to as "wiggles" ) of the perfect brane worldsheet 
will be given, according to the characteristic equation (3.10), by 

c E 2 = | ■ (4-10) 

By substituting (4.8) in the (3.3), and using the expression (1.16) for the acceleration w M of the 
preferred unit vector (4.8), the extrisic equation governing the free motion of the support surface 
of any (proper) perfect brane is reducible to the standard formfl[p4[ 

c E 2 K» = (l-cl)±H u u» . (4.11) 

In the Dirac-Goto-Nambu case one has U = T and hence c E 2 = 1, (i.e. the "wiggle" speed is that of 
light) so that the right hand side of (4.11) will drop out, leaving the harmonicity condition that is 
expressed by the vanishing of the curvature vector given by (1.19). At the opposite extreme from 
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this relativistic limit case, one has applications to such mundane examples as that of an ordinary 
soap bubble membrane or violin string for which one has c E 2 << 1. 

Except in the Dirac-Goto-Nambu limit case, the extrinsic equation of motion (4.12) will need 
to be supplemented by the dynamic equations governing the evolution of the internal fields, and in 
particular of T, U, and u^, on the brane surface. The simplest non trivial possibility is the case of 
a perfect brane model that is "barotropic" , meaning that its tension T (or equivalently, in the case 
of a fluid, its pressure P) is a function only of the energy density U in the preferred rest frame. 
In this barotropic or "perfectly elastic" case, the complete set of dynamical equations governing 
the internal evolution of the perfect brane will be provided just by the tangential projection of the 
force balance law (2.24) which, in the free case to which the discussion of the present section is 
restricted, has the simple form 

V X v %T fUJ = . (4.12) 

In the barotropic case (which includes the kind of cosmic string models |3lJ j3^] [|?4j appropriate 
for describing the Witten type superconducting vaccuum vortices whose investigation provided 
the original motivation for developing the kind of analysis presented here) experience of the "im- 
proper" case of an (n-l)-brane, i.e. the familiar example of an ordinary barotropic perfect fluid, 
suggests the convenience of introducing an idealised particle number density, v say, and a cor- 
responding (relativistic) chemical potential, or effective mass per idealised particle, /i say, that 
are specified modulo an arbitrary normalisation factor by the equation of state according to a 
prescription of the form 

lnu= , In a = / , (4.13) 

J U-T ' ' J T-U' 1 ' 

which fixes them modulo a pair of constants of integration which are not allowed to remain inde- 
pendent but are related by the imposition of the restraint condition 

fj,u = U — T . (4.14) 

The quantities introduced in this way can be used for defining a surface particle current density 
vector v p and a corresponding dynamically conjugate energy - momentum (per particle) covector 
Hp that are given by 

v p = vu p , Hp = fiUp , (4-15) 

in terms of which the generic expression (4.9) of the surface stress momentum energy density tensor 
of a perfect brane can be rewritten in the form 

T" a = vP^a - Tr,P a , (4.16) 

whose divergence can be seen to be given identically by 

%TP a = h*%v p + 2^"V [pMff] - TK a . (4.17) 
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In the free case for which this divergence vanishes, projection orthogonal to the worldsheet gives 
back the extrinsic dynamical equation (4.11), while by contraction with v p one obtains the simple 
surface current conservation law 

%v p = Q . (4.18) 

The remainder of the system of equations of motion of the perfect brane can hence be obtained 
from (4.12), as a surface generalised version of the standard perfect fluid momentum transport 
law|5^|, in the form 

<XV[^ P ] = , (4.19) 

which is the integrability condition for the existence of a surface potential function, ip say, such 
that /i p = W p (p. 

Appart from the extrinsic perturbations of the world sheet location itself, which propagate with 
the "wiggle" speed c E (relative to the frame deterined by u M ) as already discused, the only other 
kind of perturbation mode that can occur in a barytropic string are longitudinal modes specified 
by the variation of U or equivalently of T within the world sheet. As in the "improper" special case 
of an ordinary perfect fluid, for which they are interpretable just as ordinary sound waves, such 
longitudinal "woggle" perturbations can easily be seen from (4.16) and (4.17) to have a relative 
propagation velocity c L say that is given^]|34 by 



2 vdfx dT 



In order for a barotropic model to be well behaved, c E 2 and c L 2 must of course both be positive, in 
order for the velocities to be real (as a condition for local stability), and they must also both be 
less than unity in order to avoid superluminal propagation (as a condition for local causality) . An 
important qualitative question that arises for any proper barotropic p-brane, i.e. one with 1< p<n- 
1 (excluding the trivial case p=0 of a point particle, and at the other extreme the "improper" case 
p=n-l of a fluid) is whether the extrinsic propagation speed is subsonic, meaning that "wiggles" 
go slower than "woggles" , i.e c E < c L , transonic, meaning that the speeds are the same, c E = c L 
or supersonic, meaning that "wiggles" go faster than "woggles", c E > c L . Much of the earliest 
workf39[ jGOf fil[ f02| [gp Q [[sH] on the Witten type superconducting vaccuum vortex phenomena 
was implicitly based on use of a string model of subsonic type, but a more careful analysis of 
the internal structure of the Witten vortex by Peter Jg(| (g^| |i8j has since shown that models of 
supersonic type are more appropriate. On the other hand such everyday examples as that of an 
ordinary violin string are subsonic, a condition that is sufficient fl6l| though not necessary for 
stability of the corresponding circular centrifugally supported loop configurations f7lf |72]| [f73f |74]| 
that will be discussed later on. 

A particularly important special case is that of a barytropic p-brane of permanently transonic 
type, meaning one characterised by an equation of state of the non dispersive constant product 
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form 

UT = m 2p+2 , (4.21) 

where to is a constant having the dimensions of a mass, for which application of the formulae 
(4.10) and (4.18) gives c E 2 = c L 2 not just for some critical transition value of U or T (as can occur 
for other kinds of equation of state) but for all values. In the string case, p=l, it can easily be 
demonstrated that this non dispersive permanently transonic model is represents the outcome 
of a (rather artificial) dimensional reduction process first suggested by Nielsen |7a| J76| [7^] ||78[ . As 
will be described below, this model can be shown|79| to be governed by equations of motion that 
are explicitly integrable in a flat empty background, and it can be used]79| p0[ to provide what 
(contrary to a misleading claim|8lj that has been made) can be a highly accurate description of 
the macroscopic effect of "wiggles" in an underlying cosmic string of the simple Goto-Nambu type. 



5 Essentials of Elastic String Dynamics 

Between the hypersurface supported case of a membrane and the curve supported case of a point 
particle, the only intermediate kind of brane that can exist in 4-dimensions is that of 1-brane, i.e. 
a string model, which (for any background dimension n) will have a first fundamental tensor that 
is expressible as the square of the antisymmetric unit surface element tensor = — given by 
(1.28): 

rf v = £» p £» v . (5.1) 

In the case of string (to which the remainder of this article will be devoted) the symmetric surface 
stress momentum energy tensor T^ v that is well defined for any "pure" (i.e. unpolarised) model 
will not only be expressible in the generic form (4.6) that expresses the (in the string case trivial) 
property of spacial anisotropy with respect to the preferred timelike or null eigenvector /3 M , but it 
will also be expressible generically in bicharacteristic form as 

= , (5.2) 

in terms of a pair of timelike or null tangent vectors f}£ that can be seen to be extrinsic bicharacteric 
vectors, meaning that they lie respectively allong the two ("right moving" and "left moving") 
directions of propagation of extrinsic perturbations, and are thus respectively orthogonal to the 
covectors x M given by the extrinsic characteristic equation (3.10) whose solutions can be seen from 
(5.2) to be given by (3^Xn = or fi^Xv = 0- The expression (5.2) does not completely determine the 
bicharacteristic vectors /3£, but leaves open the possibility of a reciprocal multiplicative rescaling 
of their magnitudes, subject to the condition that their scalar product remains invariant, its value 
being obtainable by comparison of (5.2) with (4.6) as 

P + ^ = T V V =-{U + T) . (5.3) 
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The geometric mean of the magnitudes of the bicharacteristic vectors remains similarly invariant, 
with value given by the magnitude of the preferred eigenvector (3^ as given by (4.7), i.e. 

(/5+p/VO {P-vP-) = {W? = ( U ~ T ? ■ ( 5 - 4 ) 

To obtain the equations of motion of the string we need to use the force balance equation (2.25), 
which just takes the form 

%T^=T, (5.5) 

for a string that is isolated, where ff (which simply vanishes in the free case) is the force exerted by 
any background fields, such as the electromagnetic field and ambient (Higgs or other) fluid that are 
allowed for in the explicit expression (2.26). (The only difference if the string were not isolated, 
i.e. if it belonged to the boundary of one or more attached membranes, is that an additional 
contact force contribution ft 1 of the form (2.25) would also be needed as well as the background 
contribution ft 1 .) Using the Weingarten integrability condition (1.13) which is equivalent in the 
string case to the the projected Lie commutativity condition 

1F V (/V%9_" - P^%P + V ) = , (5.6) 

it can be seen from (5.2) that the extrinsic equations of motion governing the evolution of the 
string worldsheet will therefore [Q be expressible in characteristic form as 

1P V r , (5.7) 

for either choice of the sign ±, one version being obtainable from the other and vice versa by (5.6). 

Except where one of the bicharacteristic vectors (3 J 1 is null, which by (5.5) will only occur 
where U — T (as is everywhere the case for the Goto-Nambu model, for which both bicharacteristic 
vectors are null everywhere) it will be possible to fix their normalisation by requiring that they 
have the same magnitude, /3 ±fl f3 ± M = — (T — U) and hence to define a corresponding pair of unit 
bicharacteristic vectors 

u ± f = (U-T)- 1/2 [3 ± t 1 , u ± „« ± "=-l, (5.8) 

in terms of which the unit timelike eigenvector u M introduced in (4.8), together with an orthogonal 
spacelike unit eigenvalue w M , will be given by 

u»= ^"^ -(m^+m,") , v» = v/1 ~ Ce -( U+ ^- U _^) . (5.9) 

2 2c E 

Such mutually orthogonal unit eigenvectors constitute a preferred diad, in terms of which 
the antisymmetric unit surface element and the (first) fundamental tensor ■qt 11 ' of the string 
worldsheet will be given by 

^ = uV-«V, f = -«V+«V, (5.10) 
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and in terms of which the standard form (4.9) for the surface stress momentum energy density 
tensor will reduce to the simple form 

= Uu"u v - T// . (5.11) 

Proceeding directly from this last form, the extrinsic dynamical equation (5.6), i.e. the surface 
orthogonal projection of (5.5), can be rewritten as 

i_t (Uu v - Tv w ) =l_t T , (5.13) 

using the notation 

ii p = u v V v u v , = v v V„v v . (5.14) 
in terms of which the worldsheet curvature vector is expressible as 

K» =±t (v' v - ii v ) . (5.15) 

These last formulae illustrate a special feature distinguishing string models from point particle 
models on one hand and from higher dimensional brane models on the other, namely the dual 
symmetry ]3l[ that exists at a formal level between the timelike eigenvector u p and the associated 
eigenvalue U on one hand, and the the dynamically conjugate quantities that are the spacelike 
eigenvector and v p and the associated eigenvalue T . This feature is particularly informative in the 
barotropic or "perfectly elastic" case discussed at the end of the previous section, i.e. in the case 
when the internal dynamics is controlled just by an equation of state giving the energy density U 
as a function of the tension T only, since in there will be a corresponding duality relation between 
the number density v and the associated chemical potential \i as introduced in (4.13), which is 
equivalent to the mutually dual pair of differential defining relations 

dll = /idv , dT = —vd\x . (5.16) 

Defining the surface duals of the current vector v p and momentum covector fi p given by (4.15) as 

-kv p — E pa v a = vv p , *fi p = £ pa n<j = ^v p , (5-17) 

one can use (4.17) to obtain the force balance equation for a barotropic string model in the self 
dual form 

t i p V a v a + *v p V a *m ct + ±- a p (Uu a - Tv" 7 ) = J p , (5.18) 

whose orthogonally projected part gives back the self dual extrinsic force balance equation (5.13), 
while its tangentially projected part gives a mutually dual pair of surface current source equations 
of the form 

(U - T)%v° = -v»~f p , (U - T)% = V7 P ■ (5.19) 



31 



It can be seen from this that, like the surface number current density v" , the dual current will 
also be conserved in a free barotropic string, i.e. when the external force f p is absent. 



It is useful) 31 to introduce a dimensionless state function called the "characteristic potential" 



■d say that is constructed in such a way as to be self dual according to the differential relation 



M = J^. (5.20) 

V 

This function is convenient for the purpose of writing the internal equations of motion obtained by 
tangential projection of (5.5) in a characteristic form that will be the analogue of the characteristic 
version (5.7) of the extrinsic equations of motion obtained by the orthogonal projection of (5.5). 
The longitudinal bicharacteristic unit vectors £ ± M say (the analogues for the internal sound type 
waves of the extrinsic bicharacteristic unit vectors it ± M introduced above) are defineable by 

* ± '*=(l-0- 1/ V±c t t;'*) . * ± A" = -1. ( 5 - 21 ) 

where c L is the longitudinal characteristic speed as given by (4.20). In terms of these, the tangen- 
tially projected part of the force balance equation (5.5) is convertible [pl| to the form of a pair of 
scalar equations having the characteristic form 

(17 - T)£/(%d T v v %u v ) = ±£^ y T , (5.22) 

from which it is directly apparent that, as asserted in the previous section, the state function c L , as 
given by (4.20), is indeed the characteristic speed of longitudinal "woggle" propagation relative to 
the preferred frame specified by the timelike eigenvector u M . When rewritten in terms of the unit 
bicharacteristic vectors it ± M that are the extrinsic analogues of l^, the corresponding extrinsic 
equations (5.7) take the form 

(U - T) it uj>% < =±t ~f , (5.23) 

(It is to be noted that the corresponding equation (32) at the end of the original discussion |3l| of 
the characteristic forms of the equations of string motion contains a transcription error whereby a 
proportionality factor that should have been y/(U — T)T was replaced by its non-relativistic - i.e. 
low tension - limit, namely y/UT.) 

What essentially distinguishes different physical kinds of "perfectly elastic", i.e. barotropic, 
string model from one another is the form of the equation of state, which depends on the microsropic 
internal details of the vacuum defect (or other structure) that the string model is supposed to 
represent. In a cosmological context, the most important special case, which will be discussed in 
detail later on, is given by the constant product formula (4.21), but for terrestrial applications the 
most useful kind of equation of state for general purposes is the one named after Robert Hooke, who 
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is recogniseably the principal founder of string theory as a branch of physics in the modern sense 
of the word. (As for supcrstring theory, it remains to be seen whether it will become established 
as a branch of physics or just as a topic of mathematics.) Hooke's famous law to the effect that 
the tension, T is proportional to the extension, i.e. the change in 1/v in the notation used here, 
provides a very good approximation in the low tension limit that is exemplified by applications 
such as that of an ordinary violin string. This law is expressible within the present scheme as 

U(Y + T) = UqY + §T 2 , (5.24) 

where Uq is a constant interpretable as the rest frame energy in the fully relaxed (zero tension) 
limit, and Y is another constant interpretable as what is known (after another pioneer worker on 
the subject) as Young's modulus. 

A more general category of non-barotropic string models is of course required, not only in such 
terrestrial applications as (above ground) electric power transmission cables (the subject of my first 
research project as an undergraduate working working for commercial industry during a vacation), 
but also for many cosmologically relevant cases that might be envisaged (such as that of a "warm 
superconducting string" [Q) in which several independent currents are present. Nevertheless most 
of the cosmological string dynamical studies that have been carried out so far have been restricted 
(not just for simplicity, but also as a very good approximation in wide range of circumstances) 
to models of the the "perfectly elastic" barotropic type, and more specifically to the category 
(including the Hookean and Goto Nambu examples as opposite extreme special cases) of string 
models describable by a Lagrangian of the form 

£ = A + JM M + W^B^ , (5.25) 

in which, in order to satisfy the vorticity flux conservation requirement (as derived in Section 2 
from the condition of invariance with respect to gauge changes in the background Kalb-Ramond 
field B^v) one must have 

W» v = kS" u , (5.26) 

for some fixed (quantised) circulation value K say, while another fixed coupling constant, e say, 
specifies the electromagnetic current density, in the form 

1» = ec" , (5.27) 

in terms of the only independent field variable on the world sheet, which can be taken to be a 
stream function ip which enters only via the corresponding identically conserved number density 
current 

c» = E^ Pv , p v = V„V , (5.28) 

whose squared magnitude, 

X = c^cf = -p»p v = -UH^i , (5.29) 
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is the only argument in a (generically non-linear) master function. 



A = A{ x } . (5.30) 

Previous studies of cosmic strings have mostly been restricted to cases involving allowance for a 
Kalb-Ramond coupling (which is relevant for vortex defects of "global type" p^||^3|||8^] in the 
absence of an internal current, or alternatively to cases (such as are exemplified |3l|][[32) by Wit- 
ten's local superconducting string mechanismQ) involving allowance for an internal current in the 
absence of Kalb-Ramond coupling. However there is no particular difficulty in enviseaging the 
presence of both kinds of coupling simultaneously as is done here. A new result that is obtained 
thereby (in Section 6) is that for stationary and other symmetric configurations, the effect of the 
coupling to tht Kalb Ramond background is expressible as that of a fictitious extra electromagnetic 
field contribution. 

It is the specific form of the master function that determines the relevant equation of state - 
or if its range of validity is sufficiently extensive, the pair of equations of state, corresponding to 
the qualitatively different regimes distinguished by the positivity or negativity of x ~ by which the 
dynamics of the model is governed. The derivative of the master function provides a quantity 

K = 2^ , (5.31) 

dx 

in terms of which the surface stress momentum energy density tensor defined by (2.7) will be given 
by 

= /CpV + An^ . (5.32) 

This only depends on the internal field gradient and on the imbedding and the background ge- 
ometry: it involves neither the background electromagnetic and Kalb Ramond fields nor their 
respective coupling constants e and k. The role of the latter is just to determine the background 
force density, which will be given, according to the general formula (2.26), by 

7 p = eF w c" + ^N mv £» v . (5.33) 

Unlike the force density, the stress momentum energy density tensor, as given by (5.32), is formally 
invariant under a duality transformation Q [[H} whereby the variational momentum covector p^ , the 
current c M , and the master function, A itself, are interchanged with their duals, as denoted by a 
tilde, which are definable by 

c' 1 = JCp^ , p^ = /Cc M , x = c M c M = -p v p v , (5.34) 

and by the relation 

XJC = A - A = -xK , (5.35) 
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with 

t = 2^ = IC- 1 , x= -IC 2 x ■ (5.36) 

ax 

In terms of these quantities the stress momentum energy density tensor can be rewritten in either 
of the equivalent mutually dual canonical forms 

c»p v + Arf v = T"„ - (fp v + hf v . (5.37) 

Substitution of these respective expressions in (5.5) then gives the force balance equation in the 
equivalent mutually dual forms 

Pp %c p + 2£ p V [pPp] + AK P =J p = p p %cf + 2c»%p p] + AK P . (5.38) 

The independent internal equations of motion of the system are obtainable by contraction with the 
independent mutually orthogonal tangent vectors p p and p p , the latter giving what is just a kine- 
matic identity, since, by construction according to the prescription given above, one automatically 
has 

V^ 1 = = £ WV^Pp , p p f p = . (5.39) 

The only part of the internal equations of motion that is properly dynamical (from the point of 
view of the variation principle with respect to A) is thus just the part got by contraction with p p , 
which gives 

- X %& = P P 1 P = x£ pfl %P P ] , P P 7 P = e X £ pp F pp . (5.40) 

It is to be noted that unlike the electromagnetic force, the Joukowski force (interpretable as a 
manifestation of the Magnus effect) arising from the Kalb-Ramond coupling in (5.33) will always 
act orthogonally to the world sheet and thus will not affect the internal dynamics of the string. It 
is evident that the internal force balance equation (5.40) can be rewritten simply as 

£ pp (2V [p p p] +eF pp ) = , (5.41) 

which can be construed as the (Poincare type) integrability condition for the local existence of 
a scalar tp on the worldsheet in terms of which the dual momentum will simply be given as the 
tangentially projected gauge covariant derivative: 

Pp = D pV , D p = V p p Dp , D p — V p eA p . (5.42) 

In applications to Witten type superconducting strings, this field tp will be interpretable |3l| as 
being proportional to - and therefore with respect to suitably scaled units equal to - the phase 
of an underlying complex bosonic field. The scalar tp can be used as the primary independent 
variable, instead of ip, in an alternative (dynamically conjugate) variational formulation using a 
Lagrangian constructed from A instead of A. 
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If its range of definition is sufficiently extended, the master function A will determine not just 
one but a pair of distinct equations of state, one applying to the "magnetic" regime where the 
current c p is spacelike so that \ is positive, and the other applying to the "electric" regime where 
the current c p is timclike so that the \ 1S negative. In either of these distinct regimes - though not 
on the critical intermediate "null state" locus where the current is lightlike so that \ vanishes - the 
mutually dual canonical forms (5.37) of the stress momentum energy density tensor are replaceable 
by the equivalent more manifestly symmetric standard form 

T" u = - c p c u + 4 ~c p ~c v . (5.43) 
X X 

By comparison with the corresponding expression (5.11) it can thus be seen that the physical 
interpretation - as energy density U, tension T, particle number density v and effective mass per 
particle /i - of the state functions introduced in the variational function theat has just been given 
will be given, in the "magnetic" regime where c p is spacelike (the first possibility to be considered 
in the early work on the Witten mechanismQ), by 

A = -T , A = —U , X = , X = -v 2 , £=-, (5.44) 

A 4 

with 

c p = fiv p , c p = vu p , p p — [iu p , p p — Wp . (5.45) 

On the other hand in the "electric" regime where c M is spacelike, one obtains what in general will 
be a different equation of state given parametrically by the dual set of relations 

A = —U , A = — T x = -v 2 , X = v\ £ = ( 5 ' 46 ) 

with 

c p = vu p , c p = fiv p , p p — vv p , p p — fiUp . (5-47) 

In realistic conducting string models one can not expect the appropriate master function A{%} 
to be given exactly by a simple analytic formula, but only by the output of a detailed numerical 
computation of the internal structure of the underlying vacuum defect [85| |6(| pTIj] |]68| . There 
is however an artificial mechanism (proceeding by a Kaluza Klein type projection of a simple 
Goto-Nambu model in an extended background with an extra space dimension) first suggested by 
Nielscn|75| that can easily be shownQ to provide a particularly elegant illustration of the above 
formalism with a master function given (for an appropriate choice of the scaling of the stream 
function ip) in terms of a constant mass parameter m by 

A = \J to 4 — m 2 x , A = \J to 4 — m 2 x ■ (5.48) 

This has the rather exceptional but by no means unique property of providing an equation of state, 
for the relationship between U and T, that has the same form in the magnetic regime x > > X 
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as in the electric regime x > > X- What is quite unique about this equation of state is its non- 
dispersive - i.e. permanently transonic - "constant product" form (4.21), which can be seen|79fl 
to be expressible parametrically in terms of a dimensionless self dual state function ■& of the kind 
introduced in (5.20) by 



U = m z \l 1 + -^r = m 2 coth $ , T = m 2 Jl-4 = ?7i 2 tanh $ , (5.49) 

V m 2 V m 2 

which gives 

c E = c L = tanh$ , (5.50) 

so that by (4.21) and (4.20), the extrinsic ("wiggle") and longitudinal ("woggle") bicharacteristic 
vectors will coincide, having the form 

tj l = uj 1 = coskdu 11 ±smkdv» . (5.51) 

The revelation^] of this unique transonicity property invalidates the too hasty claim|7(| that 
Nielsen's elegant artifice|7f| effectively represents the outcome in the "pure" string limit (in which 
effects of finite vortex thickness are neglected) of the Witten mechanism Q. The latter requires a 
model Q of generically dispersive typeQ that has been shown by the work of PeterQ^Q 
to be characterised typically by supersonicity, c E > c L . The elegant transonic model can neverthe- 
less provide a crude but mathematically convenient approximation to a realistic description of a 
superconducting string that should be a considerable improvement not only on the use|8(| for this 
purpose of an unmodified Goto-Nambu model, but also on the more commonly used description 
provided by the naively linearised model|^9)|^][^l|||6^][^3|||6^][^5) with equation of state given by 
constancy of the trace (5.3) which, by (4.10) and (4.20), is evidently characterised by permanent 
subsonicity, c E < c L = 1. 

In terms of the bicharacteristic vectors (5.51), the intrinsic equations of motion (5.19) can be 
recombined|13| [p2| as an equivalent pair of divergence relations in the form 



v P ((c/ 



T)£ ± p )=-L p f p , (5.52) 



which shows that the "left" and "right" moving "bicharacteristic currents", (U — T)£ + p and (U — 
T)£_ p will each be conserved separately in the free case, i.e. when / vanishes. As an alternative 
presentation of the tangential force balance equations for this permanently transonic model, a little 
algebra suffices to show that the intrinsic equation of motion (5.22) can be rewritten for this case, 
in a form more closely analogous to that of the corresponding extrinsic equation of motion (5.23), 
as the pair of equations 

(U - T)if v £/V p l T » = « + i^i^Yr ■ (5.53) 
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This tangentially projected part of the dynamic equations can now be recombincd [|79| with its 
orthogonally projected analogue (5.23), so as to give the complete set of force balance equations 
for the non dispersive permanently transcharacteristic string model (4.21) as the extremely useful 
pair of bicharacteristic propagation equations 

(U - T)£/V p £/ = (g^ + W)7- • ( 5 - 54 ) 

The advantage of the characteristic formulation (5.54) is that it takes a particularly simple form 
when expressed in terms of the corresponding characteristic coordinates a ± on the worldsheet as 
defined by taking the "right moving" coordinate er+ to be constant allong "left moving" charac- 
teristic curves, and taking the "left moving" coordinate a- to be constant along "right moving" 
characteristic curves, with the convention that the correspondingly parametrised bicharacteristic 
tangent vectors £^ = dx fJ, /da ± should be future directed. In the case of a free motion (i.e. when 
the force term on the right of (5.54) vanishes) in a flat background, this simplification can be 
used j79| to obtain the complete solution of the dynamical equations in a very simple explicit form. 
The way this works is that (5.54) reduces just to 

dl£/da w =0, (5.55) 

whose general solution is given in terms of a pair of generating curves x^{o~} as a sum of single 
variable functions by the ansatz 

x^ = x£{<t+} + x^{ct-} , (5.56) 

which gives t£ = x£ using the dot here to denote the ordinary derivatives of the single variable 
functions with respect to the corresponding characteristic variables. The solution (5.56) generalises 
a result that is well known for more familiar but degenerate Goto Nambu case, in which the 
tangents to the generating curves are required to be null, £+£ ±IJ , = (so that with the usual 
normalisation their space projections lie on what is known as the Kibble Turok sphere |87|| ). The 
only restriction in the non - degenerate case is that they should be non-spacelike and future directed 
(the corresponding projections thus lying anywhere in the interior, not just on the surface, of a 
Kibble Turok sphere) the unit normalisation condition (5.21) being imposable as an option, not an 
obligation, by choosing the parameter a to measure proper time allong each separate generating 
curve. 

This special property of being soluble by an ansatz of the same form (5.56) as has long been 
familiar for the Goto Nambu case can immediately be used to provide a new direct demonstration of 
the validity of the non-dispersive constant product model (4.21) for describing the average motion of 
a "wiggly" Goto Nambu string. The previously unpublished justification presented here is needed 
because my original argument |Q was merely of a qualitative heuristic nature, while Vilenkin's 
mathematical confirmation [p0[ was based on indirect energetic considerations, and has been called 
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into question |yj on the grounds that it did not cover the most general class of wiggles that can 
be envisaged. It has recently been shown by Martinj8£| that Vilenkin's method|8l| (extending to 
the tension T the concept of "renormalisation" that had previously been introduced for the energy 
density U by Allen and ShellardQ (9(J) can in fact be generalised straightforewardly so as to 
confirm the validity of (4.21) for all kinds of "wiggle" perturbations (subject only to the restriction 
that their amplitudes should not be so large as to bring about a significant rate of self intersection) . 
Though based on manifestly muddled reasonning, the purported contrary demonstration |yj of 
higher order "deviations" from the constant product form (4.21) for the effective equation of state 
has left a residue of controversy and confusion. The need for an absolutely clear refutation of such 
allegations has motivated the formulation, as a more direct alternative to the energetic analysis 
developed by Vilenkin and Martin, of the new derivation presented here. 

The new justification for the use of the elastic string model characterised by (4.21) consists 
simply of the observation that such a model implicitly underlies the diamond lattice discretisation 
that, since its original introduction by Smith and Vilenkin |9l|], has been commonly employed by 
numerical simulators |92j as a very convenient approximation scheme - of in principle unlimited 
accuracy - for the representation of a Goto-Nambu string worldsheet. As a way of replacing 
the exact continouous description by a discrete representation such as is necessary for numerical 
computation, the idea of the Smith Vilenkin method is simply to work with a pair of discrete sets 
of sampling points x£ T = x^{a T } determined by a corresponding discrete set of parameter values 
cr r on the generating curves of the exact representation (5.56). This provides a "diamond lattice" 
of sample points given (for integral values of r and s) by 

x is — x +r + x _s ! (5.57) 

that will automatically lie exactly on the "wiggly" Goto Nambu worldsheet (5.56), which is thus 
represented to any desired accuracy by choosing a sufficiently dense set of sampling parameter 
values cr r on the separate "wiggly" null generating curves x^{o~}. The new remark I wish to make 
here is simply that the chosen set of sample points x£ T — x^{<7 T } on the separate "wiggly" null 
generators can also be considered to be sample points on a pair of smoothed out, and thus no 
longer null but timelike, interpolating curves that, according to the result [79| demonstrated above, 
can be interpreted according to (5.56) as generating a corresponding solution of the equations of 
motion for an elastic string model of the kind governed by (4.21). The not so "wiggly" elastic 
string worldsheet constructed by this smoothing operation will obviously be an even better ap- 
proximation to the exact "wiggly" Goto Nambu worldheet than the original Smith Vilenkin lattice 
representation, which itself could already be made as accurate as desired by choosing a sufficiently 
high sampling resolution. No matter how far it is extrapolated to the future, the smoothed elastic 
string worldsheet generated according to (5.55) can never deviate significantly from the underlying 
"wiggly" Goto-Nambu worldsheet it is designed to represent because the exact worldsheet and the 
smoothed interpolation will always coincide precisely at each point of their shared Smith Vilenkin 
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lattice (5.57). This highly satisfactory feature of providing a potentially unlimited accuracy could 
not be improved but would only be spoiled by any "deviation" from the originally proposed Q 
form (4.21) for the effective equation of state. 

After thus conclusively establishing that the permanently transonic elastic string model char- 
acterised by the simple constant product equation of state (4.21) (without any higher order cor- 
rections) provides an optimum description of the effect of microscopic wiggles in an underlying 
Goto Nambu model so long as self intersections remain unimportant (as was assumed in all the 
discussions ^ |80fl flsil |8S|] cited above), it remains to be emphasised that the neglect of such in- 
tersections will not be justified when the effective temperature) 79 of the wiggles is too high 
(as will presumably be the casejl3| during a transient period immediately following the string - 
forming phase transition). The result of such intersections will be the formation of microscopic 
loops, of which some will subsequently be reconnected, but of which a certain fraction will escape. 
Estimation of the dissipative cooling force density that would be needed to allow for such losses 
remains a problem for future work. 



6 Symmetric Configurations including Rings and their Cos- 
mological Implications 

Whenever the background space time metric is invariant under the action of a (stationarity, ax- 
isymmetry, or other) continuous invariance group generated by a solution fc M Killing equation 
(2.31), i.e. in the notation of (2.20) 

k£g^ = , (6.1) 

then any string that is isolated (i.e. not part of the boundary of an attached membrane) will have 
a corresponding momentum current (interpretable, depending on the kind of symmetry involved, 
as representing a flux of energy, angular momentum, or whatever) given by 

= T^k u . (6.2) 

In accordance with (2.33), this will satisfy a source equation of the form 

%P» = , (6.3) 

which means that the corresponding flux would be strictly conserved when the string were not 
just isolated but free, i.e. if the background force / were zero. When the string is subject to a 
background force of the Lorentz-Joukowsky form (5.33) that arises from background electromag- 
netic and Kalb-Ramond fields, then provided these background fields are also invariant under the 
symmetry group action generated by fc M , i.e. in the notation of (2.20) 

k£A M = , k£B MV = , (6.4) 
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it can be seen to follow that, although the physically well defined surface current will no longer 
be conserved by itself, it still forms part of a gauge dependent generalisation, V 11 say, that is strictly 
conserved, 

%V» = , (6.5) 
and that is given, in terms of the gauge dependent generalisation 

= T"„ + JM„ + W w B vp (6.6) 

of the surface stress momentum energy density tensor, by 

V" = T" v k v =7^+ {ec^A v + nE w B vp ) . (6.7) 

The purpose of this section is to consider string configurations that share the background 
symmetry under consideration. A configuration that is symmetric in this sense will be characterised 
by the condition 

_Lt k" = , (6.8) 

meaning that the symmetry generator fc M is tangential to the worldsheet, and the corresponding 
Lie invariance condition on its surface stress momentum energy density tensor will have the form 

k^T up = 2T^V p k p) . (6.9) 

Under such conditions, as well as the ordinary momentum flux P p , what may be termed the adjoint 
momentum flux, 

fpM = jfe»4r„f , iT v » = £ vp £^TP a , (6.10) 

will also obey an equation of the form (6.3), i.e. 

f V^=7^, (6.11) 

and hence would also be conserved if the string were free. In the presence of an electromagnetic 
or Kalb Ramond background, it can be seen that, like the ordinary momentum flux, this adjoint 
momentum flux has a gauge dependent extension, 

tp" = k^T^ , ^T v " = EvpE^T",, , (6.12) 

that will share with the generalised momentum flux V p (that would be conserved even if the string 
did not share the symmetry of the background) the property of obeying a strict surface current 
conservation law, namely 

% f P" = 0. (6.13) 

The group invariance conditions 

krVvPp + = , WV v pp + p^ p p v = , (6.14) 
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that are the analogues of (6.9) for the separate mutually dual pair of internal momentum covectors 
p^ and p^ associated with the internal current within the string, can be rewritten, with the aid of 
the corresponding electromagnetic background invariance condition (6.4), in the form 

k p £ pv (V v p„ + & -F vp ) = £ p ^ v lo , We^VvPv = £ pv V„P , (6.15) 

where to and p are are the Bernoulli type scalars given by 

uj = p + eA^\ P=§^, P = P^- (6.16) 

If fc M is timelike, the corresponding symmetry will be interpretable as stationarity, while the more 
restrictive case[^3| in which the string is actually static (in the sense that there is no transverse 
current component relative to the background rest frame determined by k^) will be given by the 
condition that the second Bernouilli constant, (3 should vanish. 

It can be seen from (6.15) that the internal equations of motion (5.39) and (5.41) are equivalent 
in this group invariant case simply to the corresponding pair of Bernouilli type conservation laws 
to the effect that u and j3 (but, unless the electromagnetic field is absent, not (3) should both 
be constant over the worldshcet. This observation allows the problem of solving the dynamical 
for a symmetric configuration of the kind of ("perfectly elastic" , i.e. barocentric) string model 
under consideration to one of solving just the extrinsic equations governing the location of the 
worldsheet. A recent investigation^^ based on the systematic use of variational methods in the 
restricted case for which the Kalb-Ramond coupling was absent has drawn attention to the interest 
of extrapolating the Bernoulli constants outside the supporting worldsheet as a pair of scalar fields 
defined over the entire background spacetime by the uniformity conditions 

V^o; = , V p /3 = , (6.17) 

and formulating the problem in terms a certain particular worldsheet tangent vector that is defined 

by 

X n = E ^\ Pv = k v £ vp T p * . (6.18) 
This vector will be expressible in terms of the variables introduced in previous section as 

A- = ^ + &S> , (6.19) 

X X 

while the Killing vector itself will be expressible in analogous form by 

k^ = -i- + -J— c^ . (6.20) 

A-A A-A v ' 

The latter gives for the (real or pure imaginary) Killing vector amplitude, V say (which in the 
case of stationary symmetry, for which it is real, will be interpretable as an effective gravitational 
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potential field) , an expression of the form 

ft 2 ft 2 

V 2 = -kn^ = l — + 1 ^, (6.21) 

X X 

that, in view of (6.16), and of the state functional relationship between \ an d Xi can be solved 
for any particular choice of the constant "tuning parameters" u> and (3 to determine the internal 
variables \ and \ as functions of the (gravitational type) potential V and the (electric type) 
potential A^k^, and hence by implication as scalar fields over the entire background space, not 
just on the worldsheet where they were originally defined. In a similar manner, it can be seen that 
the contraction of the worldsheet generating vector X^ with the Killing vector field will be given 
simply by 

k^X„ = 0P , (6.22) 
and that its (real or pure imaginary) amplitude, X say, will given by an expression of the form 

X 2 = X^X. = P±- + P±- (6.23) 

X X 

whereby it too is implicitly defined as a function of V and A^k^, and hence as a well defined scalar 
field, not just on the worldsheet but, by (6.16), also over the background as a whole, its gradient 
being given by 

XV^X = KkVV^V + e~p(k- A)4v M (A, y fc ,y ) . (6.24) 

X 

Previous experience |Q with the case in which only the electromagnetic but not the Magnus force 
contribution is present suggests the interest of formulating the problem in terms of the propagation 
of the special generating vector X^, which (using the formulae (5.5) and (5.43) of the previous 
section) can be seen from (6.19) and (6.20) to be given by 

X v V„ X M = AA fc"^ + X v £ vp k p ~f ii , (6.25) 

where / M is the background force as given by (5.33). It can now be seen that the two preceeding 
equations can be combined to give the equation of motion for the worldsheet generating vector X* 1 
in the very elegant and convenient final form 

X V V V X» - IV M X = Tp V X v , (6.26) 

in terms of a pseudo Maxwellian field given by 

= epF^ + KN„ vp k p = 2V [v A p] , (6.27) 

where A^ is a gauge dependent pseudo-Maxwellian potential covector given by 

A„ = e/3A^ + kB^W . (6.28) 
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If fc M is timelike so that the corresponding symmetry is interpretable as stationarity, then the 
equation of motion (6.26) will be interpretable as the condition for the string to be in equilibrium 
with the given values of the constant "tuning" parameters u) and 0, (of which, as remarked above, 
the latter, j3, will vanish in the case of an equilibrium that is not just stationary but sfatacp3|]). The 
new result here is that the Joukowski type "lift" force (which was not allowed for in the previous 
analysis [Q) due to the Magnus effect on the string, as it "flies" (like an aerofoil) through the 
background medium represented by the current 3-form N pup , has just the same form as an extra 
Lorentz type electromagnetic (indeed in the stationary case purely magnetic) force contribution. 

When the genuine electromagnetic background coupling, and the similarly acting Kalb Ramond 
coupling are both absent, then (as pointed out previously |94|) the equation of motion (6.26) (that 
is the equilibrium condition for "steady flight" in the stationary case) is just a simple geodesic 
equation with respect to, not the actual background spacetime metric g M „, but the conformally 
modified metric X 2 g flu , with the conformal factor X 2 determined as a field over the background 
by (6.23) in conjunction with (6.17) and (6.21). Even when the Lorentz and Joukowski force con- 
tributions are present, the equation (6.26) governing the propagation of the world sheet generator 
X^ retains a particularly convenient Hamiltonian form, given by 

~ da ~ 97J M ' da ~ dx» ' ( ] 

for the quadratic Hamiltonian function 

H = \g^ v {H p - A^{n v - A v ) - \x 2 , (6.30) 

subject to a restraint fixing the (generically non-affine) parametrisation a of the trajectory by the 
condition that the numerical value of the Hamiltonian (which will automatically be a constant of 
the motion) should vanish, 

H = Q, (6.31) 

together with a further momentum restraint, determining the (automatically conserved) relative 
transport rate Ilji^ in accordance with the relation (6.22) by the condition 

= lo(3 . (6.32) 

The Hamiltonian momentum covector itself can be evaluated as 

77 M =X fl + A fl = £jV V , (6.33) 

or more explicitly, in terms of the original conserved generalised momentum flux as given by 
(6.7), and the gradients of the scalar (stream function and phase) potentials introduced in (5.28) 
and (5.42), as 

77 M = k v £ vp T% = V v £ Vil + luV^ + 0V^ . (6.34) 
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The advantage of a Hamiltonian formulation is that it allows the problem to be dealt with by 
obtaining the momentum covector in the form = V^S 1 from a solution of the corresponding 
Hamilton Jacobi equation, which in this case will take the form 

gT{VpS - A»)(V V S - A v ) = X 2 , (6.35) 

with X 2 given by (6.23) via (6.21), while the restraint (6.32) gives the condition 

k v V„S = (3lu . (6.36) 

Generalising results obtained previously |p3 196| for the Goto Nambu limit case, it has recently 
been shown |94j| that for a string model of the non-dispersive permanently trasonic type with the 
constant product equation of state (4.21) that is governed by the Lagrangian (5.48) (as obtained^] 
both from the Nielsen dimensional reduction mechanism and also(|79|], as explained in Section 5, 
from the more physically realistic "wiggly" string approximation) the stationary Hamilton Jacobi 
equation is exactly soluble by separation of variables in a Kerr black hole spacetime, not just of the 
ordinary asymptotically flat kind but even of the generalised asymptotically De Sitter kindj97t |98|. 
Except in the Schwarschild-De Sitter limit, where it could of course have been predicted as a 
consequence of spherical symmetry, this separability property still seems rather miraculous, re- 
flecting a "hidden symmetry" of the Kerr background that is still by no means well understood. 
The newly discovered separability property|M| is not just an automatic consequence of the sim- 
pler, though when first discovered already surprising, property of separability for the ordinary 
geodesic equation [p9[ but depends on a more restrictive requirement of the kind needed for the 
more delicate separability property of the scalar wave equation|lOO|. (It is however more robust 



than the separability properties that have turned out to hold for higher spin bosonic[101| and 
fermionic [ 102 1 [ 1Q3| [104 wave equations, and other related systems|105 [ 106 1 .) 



The most cosmologically important application of the formalism that has just been presented 
is to the equilibrium of small closed string loops in the mathematically relatively trivial case for 
which the background gravitational, electromagnetic, and Kalb-Ramond fluid are negligible, so 
that the trajectories generated by solutions of (6.26) will all just be straight lines in a Minkowski 
background. Since in such a background the Killing vector trajectories are also straight in the case 
of stationarity (though not of course for axisymmetry) this might at first be perceived implying 
that the worldsheet of a stationary string in an empty Minkowski background would necessarily be 
flat. This conclusion would exclude the possibility of closed loop equilibrium states in the absence 
of a background field, and indeed cosmologists seem (albeit for other reasons) to have entirely 
overlooked the possibility that such states might exist until the comparitively recent publication 
of an epoch making paper by Davis and Shellard[ [l07| provided the first counterexamples (the only 
previously considered equilibrium states [ 108| [6l[|[|62|]63| having been based on a magnetic support 
mechanism that was was finally judged to be too feeble to be effective except 111 as a minor 
correction). 
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The loophole in the deduction that if the trajectories generated by X^ and by fc M are both 
straight then the worldsheet must be flat is that it is implicitly based on the assumption that 
the two kinds of trajectories cross each other transversly. However there will be no restriction on 
the curvature in the transverse direction in the critical case for which the two kinds of trajectory 
coincide, i.e. for which X^ and fc M are parallel. The condition for criticality in this sense is 
expressible as 

X^E^k" = , (6.37) 

which can be seen from the original definition (6.18) to be interpretable as meaning that the tangent 
covector Xp, = £^vk v satisfies the extrinsic characteristic equation (3.10), or in other words that 
the killing vector fc M itself is bicharacteristic, in the sense of being directed allong the propagation 
direction of extrinsic perturbations of the world sheet. The criticality condition (6.37) is thus 
interpretable as a condition of characteristic flow. It means that the "running velocity" , v say, of 
relative motion of the intrinsicly preferred rest frame of the string (as determined by whichever of 
&" and B 1 is timelike) relative to the background frame specified by the (in this case necessarily 
timclikc) Killing vector fc M is the same as the extrinsic propagation velocity c E given by (4.10). 

In the presence of generic gravitational, electromagnetic and Kalb Ramond forces, the criticality 
condition (6.37) can be satisfied at particular positions, such as where there is a transition from 
a subcharacteristic running velocity, v < c E to a supercharacteristic running velocity v > c E (as 
will occur for instance on a string in a steady state of radial flow into a Schwarzschild black hole) 
or where there is a cusp, with v — c E but with subcharacteristic flow v < c E on both sides. 
However in view of the constancy of the Bernoulli "tuning" parameters ui and (3, the absence of 
any background field will always allow, and generically (the special integrable case (5.49) being 
an exception) will ensure, uniformity of the state of the string, so that the transcharacteristic flow 
condition (6.37) can be satisfied throughout its length. The space configuration of such a uniformly 
transcharacteristic steady string state can have arbitrarily variable curvature, and so is compatible 
with a closed loop topology. 

The question of closed loop equilibrium states did not arise in the earliest studies[110| of cosmic 
strings, which were restricted to the Goto-Nambu model whose bicharacteristics are always null and 
so can never be aligned with the timclikc Killing vector generating a stationary symmetry. However 
in a generic string model[30] for which the 2-dimensional longitudinal Lorentz of the internal 
structure is broken by a current, whether of the neutral kind exemplified in an ordinary violin string 
or the electromagnetic kind exemplified ]85|] p6]| |6^] p8| by Witten's superconducting cosmic string 
model, the bicharacteristic directions will generically be timelike [Q (spacelike bicharacteristics 
being forbidden by the requirement of causality) so there will be no obstacle to their alignemcnt 
with a timelike Killing vector in accordance with the criticality condition (6.37), i.e. to having 
a running velocity given by v = c E . The simple "toy" complex scalar field model on which the 
pionneering cosmic string studies fj were based had longitudinally Lorentz invariant vortex defects 
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(of "local" or "global" type depending on the presence or absence of coupling to a gauge field) 
that were describable at a microscopic level by string models (with Kalb Ramond coupling in 
the "global" case) that were indeed of the special Goto-Nambu type. However the extra degrees 
of field freedom (starting with the additional scalar field introduced by Witten in his original 
superconducting example 0) that are needed in successively more realistic models) 111 [ 1 12 1 make it 
increasingly difficult to avoid the formation of internal structure breaking the longitudinal Lorentz 
invariance and reducing the extrinsic characteristic velocity to the subluminal range c E < 1 at 
which stationary equilibrium with the critical running velocity v = c E < 1 becomes possible. 

The cosmological significance of this is that whereas Goto Nambu string loops cannot ultimately 
avoid gravitationally or otherwise radiating away all their energy 

mm mil ) 

since they have 

no equilibrium states into which they might settle down, on the other hand more general kinds 
of strings, whose occurrence would now seem at least as plausible) 111 [112|, can leave a relic 
distribution of stationary loop configurations that may survive indefinitely [ 107 1 ( 116 1 [ 117 ] [ 118 1 . One 
would expect such configurations to be those that minimise the energy for given values of the 
relevant globally conserved quantities, of which there might be a considerable number in the more 
complicated multiply conducting models that might be considered, but of which there are only a 
single pair in the "barotropic" type string considered here, namely the stream function winding 
number <f dip and the phase winding number <f> dip whose respective constancy results from the 
conservation of the mutually dual pair of currents c M and c M . 

The potential cosmological importance of such a distribution of relic loops was first pointed 
out by Davis and Shcllard| 107| |[116|[117], who emphasised that in the case of the "heavy" cosmic 
strings whose existence had been postulated to account for galaxy formation, the ensuing relics, 
even if formed with very low efficiency, would be more that sufficient to give rise to a catastrophic 
cosmological mass excess of the kind first envisaged as arising from the formation of monopolcs. 
According to a more recent and detailed order of magnitude estimate of my own[ 118 |, the di- 
mensionless coupling constant Gm 2 ss 10~ 6 characterising "heavyweight" strings, meaning those 
arising from G.U.T. symmetry breaking (for which the relevant Higgs mass scale m is within a 
factor of order a thousand of the Plank mass) would have to be reduced below a value given very 
roughly by Gm 2 « 10~ 26 to avoid a cosmological mass excess today. This estimate |l!8| (based 
on considerations of the kind discussed during the present meeting by Zurek) should be regarded 
as provisional, pending the more deeper investigation that would seem to be needed. What is 
remarkable about this tentative limit is that it corresponds to a Higgs mass scale m of roughly the 
same order as that at which electroweak symmetry breaking is believed to occur, which is rather 
suggestive in view of the fact that although the "standard" electroweak model does not give rise 
to stable string like vortex defects[119], nevertheless such defects do occur in many of its most 



commonly considered competitors[lll 



The cosmological implications (as an argument against "heavy" cosmic string formation, or 
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more positively as source of "dark matter" in the form of "lightweight" string loops due to elec- 
troweak symmetry breaking) of the long term survival of cosmic string loops, motivates more 
thorough investigation of equilibrium states that may be involved, a particularly important ques- 
tion being that of their stability. Prior to the derivation of the general symmetric string generator 
equation (6.26), the only closed loop equilibrium states to have been considered were the circular 
kind referred to by Davis and Shellard[ 107| [11€ ] [ 117 as "vortons" a term that is more appropri- 
ate for axisymmetric ring states in the "global" case than in the gauge coupled "local" case for 
which the string description is most accurate. The first general investigation^ [ 120 1 of such cir- 
cular ring states showed that under conditions of purely centrifugal support (neglecting possible 
electromagnetic corrections of the kind evaluated more recently| |l09| ) the condition (6.37) for equi- 
librium, namely the requirement of a transcharacteristic rotation speed v = c E , is such that the 
ring energy is minimised with respect to perturbations preserving the circular symmetry. However 
a more recent investigation of non-axisymmetric perturbations has shown that although there are 
no unstable modes for states of subsonic rotation]^ (as exemplified by a cowboy's lassoe loop) 
with v = c E < c L , instability can nevertheless occur for rotation in the supersonic regime [|70[ that 
(contrary to what was implicitly assumed in earlier work) 59 [ 30| |6l| |32| |6j| |34j |65| using the sub- 
sonic type of model given by a linear equation of state for which the sum U + T is constant) has 
been shown by Peter|Q|6^]|^8| to be relevant in the kind of cosmic vortex defects that have been 
considered so far. Although the first category of string loop equilibrium states to have been studied 
systematically has been that of circular ring configurations [120 , it has been made clear by recent 
workp4| ]69]| that, as explained above, arbitrary non circular equilibrium states are also possible. 
The stability such more general equilibrium states has not yet been investigated, carried out, but 
it seems plausible that while some of them may be destroyed by the recently discovered classi- 
cal instability mechanism []69| ]70[| , and also the kind of quantum tunnelling instability mechanism 
considered by Davis [1 16 1, it does not seem likely that such mechanisms could be so consistently 
efficient as to prevent the long term survival of a lot of other stationary loop states. 



Acknowledgements. 

I wish to thank B. Allen, C. Barrabes, U. Ben-Ya'acov, A-C. Davis, R. Davis, V. Frolov, 
G. Gibbons, R. Gregory, T. Kibble, K. Maeda, X. Martin, P. Peter, T. Piran, D. Polarski, M. 
Sakellariadou, P. Shcllard, P. Townsend, N. Turok, T. Vachaspati, and A. Vilenkin, for many 
stimulating or clarifying discussions. 

References 

[1] P.A.M. Dirac, Proc. Roy. Soc. London A268, 57 (1962). 



48 



[2] P.S. Howe, R.W. Tucker, J. Phys. AlO, L155 (1977). 

[3] A. Achucarro, J. Evans, P.K. Townscnd, D.L. Wiltshire, Phys. Lett. 198 
B, 441 (1987). 

[4] B. Carter, in Formation and Evolution of Cosmic Strings, ed. G. Gibbons, 
S. Hawking, T. Vachaspati, ppl43-178 (Cambridge U.P., 1990). 

[5] T.W.B. Kibble, J. Phys. A9, 1387 (1976). 

[6] A. Vilenkin, A.E. Everett, Phys. Rev. Lett 48, 1867 (1982). 

[7] E. Witten, Nucl. Phys. B249, 557 (1985). 

[8] B. Carter, J. Geom. Phys. 8, 53 (1992). 

[9] B. Carter, J. Class. Quantum Grav. 9, 19 (1992). 

[10] L.P. Eisenhart, Riemannian Geometry (Princeton U.P., 1926, reprinted 
1960). 

[11] J. Stachcl, Phys. Rev. D21, 2171 (1980). 

[12] J. A. Schouten, Ricci Calculus (Springer, Heidelberg, 1954). 

[13] B. Carter, M. Sakellariadou, X. Martin, Phys. Rev. D50, 682 (1994) 

[14] A. Vilenkin, Phys. Rev. D43, 1060 (1991). 

[15] R. Penrose, W. Rindler, Spinors and Space-Time (Cambridge U.P., 1984). 

[16] I. Bars, C.N. Pope, Class. Quantum Grav., 5, 1157 (1988). 

[17] P. Sikivic, Phys. Rev. Lett. 48, 1156 (1982). 

[18] A. Vilenkin, A.E. Everett, Phys. Rev. Lett. 48, 1867 (1982). 

[19] E.P.S. Shellard, in Formation and Evolution of Cosmic Strings, cd. G. 
Gibbons, S. Hawking, T. Vachaspati, ppl07-115 (Cambridge U.P., 1990). 

[20] Y. Nambu, Nucl. Phys. B130, 505 (1977). 

[21] N.S. Manton Phys. Rev. D28, 2019 (1983). 

[22] T. Vachaspati, A. Achucarro, Phys. Rev. D44, 3067 (1991). 

[23] T. Vachaspati, M. Barriola, Phys. Rev. Letters 69, 1867 (1992). 

49 



[24] A. Dabholkar, J.M. Quashnock, Nucl. Phys. B333, 815 (1990). 

[25] R.A. Battye, E.P.S. Shellard, "Global string radiation", DAMTP preprint 
(Cambridge, 1994). 

[26] B. Carter, Class. Quantum Grav. 11, 2013 (1994). 

[27] A. Vilenkin, T. Vachaspati, Phys. Rev. D35, 1138 (1987). 

[28] R.L. Davis, E.P.S. Shellard, Phys. Rev. Lett. 63, 2029 (1989). 

[29] M. Sakellariadou, Phys. Rev. D44, 3767 (1991). 

[30] U. Ben-Ya'acov, Nucl. Phys. B382, 597 (1992). 

[31] B. Carter, Phys. Lett. B224, 61 (1989). 

[32] A.L. Larsen, Class. Quantum Grav. 10L, 35 (1993). 

[33] J. Garriga, M. Sakellariadou, Phys. Rev. 48, 2502 (1993). 

[34] B. Carter, Phys. Lett. B228, 466 (1989). 

[35] B. Carter, Phys. Rev. 48, 4835 (1993) 

[36] J. Garriga, A. Vilenkin, Phys. Rev. D47, 3265 (1993). 

[37] A.L. Larsen, V. Frolov, Nucl. Phys. B414, 129 (1994). 

[38] J. Garriga, A. Vilenkin, Phys. Rev. D44, 1007 (1991). 

[39] J. Garriga, A. Vilenkin, Phys. Rev. D45, 3469 (1992). 

[40] J. Guvcn, Phys. Rev. D48, 4464 (1993). 

[41] J. Guvcn, Phys. Rev. D48, 5563 (1993). 

[42] B. Carter, Proc. Roy. Soc. Lond. A372, 169 (1980). 

[43] A. Polyakov, Nucl. Phys. B268, 406 (1986). 

[44] K.I. Maeda, N. Turok, Phys. Lett. B202, 376 (1988). 

[45] R. Gregory, Phys. Lett. B206, 199 (1988). 

[46] R. Gregory, Phys. Rev. D43, 520 (1993). 

[47] D. Garfinkle, R. Gregory, Phys. Rev. D41 1889 (1990). 

50 



[48] R. Gregory, D. Haws, D. Garfinklc, Phys. Rev. D42 343 (1991). 

[49] V. Silveira, M.D. Maia, Phys. Lett A174, 280 (1993). 

[50] C. Barrabes, B. Boisscau, M. Sakellariadou, Phys. Rev. D49, 2734 (1994). 

[51] B. Carter, R. Gregory, Curvature corrections to dynamics of domain walls, 
DAMTP preprint (Cambridge, 1994). 

[52] P.S. Letelier, Phys. Rev. D41, 1333 (1990). 

[53] D. H. Hartley, R.W. Tucker, in Geometry of Low Dimensional Manifolds, 
1 ( L.M.S. Lecture Note Series 150, ed. S. Donaldson, C. Thomas (Cam- 
bridge U.P., 1990). 

[54] H. Arodz, A. Sitarz, P. Wcgrzyn, Act. Phys. Polon. B, 22, 495 (1991); 23, 
53 (1992). 

[55] B. Boisseau, P.S. Letelier, Phys. Rev. D46, 1721 (1992). 

[56] B. Carter, Equations of motion of a stiff geodynamic string or higher 
brane, Meudon preprint (to appear in Class. Quantum Gravity 11, 1994). 

[57] S.W. Hawking, G.F.R. Ellis, The Large Scale Structure of Spacetime 
(Cambridge U.P., 1973). 

[58] B. Carter, in Active Galactic Nuclei, eds. C. Hazard, S.Mitton, pp 273-300 
(Cambridge U. P., 1976). 

[59] D.N. Spergel, T. Piran, J. Goodman, Nucl. Phys. B291, 847 (1987). 

[60] A. Vilcnkin, T. Vachaspati, Phys. Rev. Lett. 58, 1041 (1987) 

[61] E. Copeland, M. Hindmarsh, N. Turok, Phys. Rev. Lett. 58, 1910 (1987) 

[62] E. Copeland, D. Haws, M. Hindmarsh, N. Turok, Nucl. Phys. B306, 908 
(1988). 

[63] D. Haws, M. Hindmarsh, N. Turok, Phys. Lett. B209, 225 (1988). 
[64] D.N. Spergel, W.H. Press, R.J. Scherrer, Phys. Rev. D39, 379 (1989) 
[65] P. Amsterdamski, Phys Rev. D39, 1534 (1989). 
[66] P. Peter, Phys. Rev. D45, 1091 (1992). 



51 



[67] P. Peter, Phys. Rev. D46, 3335 (1992). 

[68] P. Peter, Phys. Rev. D47, 3169 (1993). 

[69] B. Carter, X. Martin, Ann. Phys. 227, 151 (1993). 

[70] X. Martin, Meudon preprint (to appear in Phys. Rev. D, 1994). 

[71] R.L. Davis, E.P.S. Shellard, Phys. Lett B209, 485 (1988). 

[72] R.L. Davis, E.P.S. Shellard, Nucl. Phys. B323 209 (1989) 

[73] B. Carter, Phys. Lett. B238, 166 (1990). 

[74] P. Peter, Phys. Lett. B298, 60 (1993). 

[75] N.K. Nielsen, Nucl. Phys. B167, 248 (1980). 

[76] N.K. Nielsen, P. Olcscn, Nucl. Phys. B291, 829 (1987). 

[77] A. Davidson, A.K. Wali, Phys. Lett. B213, 419 (1988). 

[78] A. Davidson, A.K. Wali, Phys. Rev. Lett. 61, 1450 (1988). 

[79] B. Carter, Phys. Rev. D41, 3869 (1990). 

[80] A. Vilenkin, Phys. Rev. D41, 3038 (1990). 

[81] J.Hong, J.Kim, P. Sikivic, Phys. Rev. Lett. 69, 2611 (1980). 

[82] B. Carter, Nucl. Phys. B412, 345 (1994). 

[83] E. Witten, Phys. Lett. B158, 243 (1985). 

[84] R.L. Davis, E.P.S. Shellard, Phys. Lett. B214, 219 (1988). 

[85] A. Babul, T. Piran, D.N. Spergel, Phys. Lett. B202, 207 (1988). 

[86] D.N. Spergel, W.H. Press, R.J. Scherrer, Phys. Rev. D39, 379 (1989). 

[87] T.W.B. Kibble, N. Turok, Phys. Lett. B116, 141 (1982). 

[88] X. Martin, Meudon preprint (1994). 

[89] B. Allen, E.P.S. Shellard, Phys. Rev. Lett. 64, 119 (1990). 



52 



[90] E.P.S. Shcllard, B. Allen, in Formation and Evolution of Cosmic Strings, 
cd. G. Gibbons, S. Hawking, T. Vachaspati, pp421-448 (Cambridge U.P., 
1990). 

[91] A.G. Smith, A. Vilenkin, Phys. Rev. D36, 990 (1987). 

[92] A. Albrecht, in Formation and Evolution of Cosmic Strings, ed. G. Gib- 
bons, S. Hawking, T. Vachaspati, pp403-419 (Cambridge U.P., 1990). 

[93] B. Carter, Class, and Quantum Grav. 7, L69 (1990). 

[94] B. Carter, V.P. Frolov, O. Heinrich, Class, and Quantum Grav. 8, 135 
(1991). 

[95] V.P. Frolov, V.D. Skarzhinsy, A.I. Zelnikov, O. Heinrich, Phys. Lett. 224, 
225 (1989). 

[96] B. Carter, V.P. Frolov, Class. Quantum Grav. 6, 569 (1989). 

[97] B. Carter, in Black Holes (Les Houches, 1972), cd. C. and B.S. DeWitt, 
pp 57-214 (Gordon and Breach, New York, 1973). 

[98] G.W. Gibbons, S.W. Hawking, Phys. Rev. D15, 2738 (1976). 

[99] B. Carter, Phys. Rev. 174, 1559 (1968). 
[100] B. Carter, Commun. Math. Phys. 99, 563 (1968). 
[101] S.A. Tcukolsky, Astroph. J. 185, 283 (1973). 
[102] W. Unruh, Phys. Rev. Lett. 31, 1265 (1973). 
[103] S. Chandrasekhar, Proc. Roy. Soc. Lond. A349, 571 (1976). 
[104] R. Guven, Phys. Rev. D22, 2327 (1980). 
[105] J-A. Marck, Proc. Roy. Soc. Lond. A385, 431 (1983). 
[106] B. Carter, J. Math. Phys. 28, 1535 (1987). 
[107] R.L. Davis, E.P.S. Shellard, Phys. Lett. B209, 485 (1988). 
[108] J. Ostriker, C. Thompson, E. Witten, Phys. Lett B180, 231 (1986) 
[109] P. Peter, Phys. Lett. B298, 60 (1993). 
[110] T.W.B. Kibble, Phys. Rep. 67, 183 (1980). 



53 



[Ill] P. Peter, Phys. Rev., D46, 3322 (1992). 

[112] W. B. Perkins, A-C. Davis, Nucl. Phys. B406, 377 (1993). 

[113] N. Turok, Nucl. Phys. B242, 520 (1984). 

[114] T. Vachaspati, A. Vilenkin, Phys Rev. D31, 3035, (1985). 

[115] R. Durrcr, Nucl. Phys. B328, 238 (1989). 

[116] R.L. Davis, Phys. Rev. D38, 3722 (1988). 

[117] R.L. Davis, E.P.S. Shcllard, Nucl. Phys. B323, 209 (1989). 

[118] B. Carter, Ann. N.Y. Acad. Sci. 647, 758 (1991). 

[119] M. James, L. Perivolaropoulos, T. Vachaspati, Nucl. Phys. B395, 534 
(1993). 

[120] B. Carter, Phys. Lett. B238, 166 (1990). 



54 



